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{57) Abstract 

Apparatus for determining a displacement characterisdc of an object comprises means (Txl, Tx2, 3, 4) for transmitting, at a given 
transmission time, a probe signal towards the object; means (5, 6) for receiving the probe signal returned by the object; means (DLl 1. 
DLl 2, PG28, PG30) for generating a detectt<»i timing signal at a delay after the transmission time, corresponding to a selected range for 
the object; first detecting means (Rx7, Rx8), coupled to the receiving means and responsive to the timing signal, for detecting the returned 
probe signal occurring at the delay of the timing signal; and second detecting means (Rxl9, Rx20) for detecting the reladve timing of the 
detection timing signal and of a signal having a predetemixned timing with respect to the transmission time; whereby a measure of the 
range of the object can be determined from the reladve timing of these two signals, from which lange measure the object displacement 
characteristic can be determined. Apparatus for transmitting and receiving data is also disclosed, as further is apparatus for assessing the 
^jproach of an object to a specified location. Analogous methods are also disclosed. 
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5 

APPARATUS AND METHOD FOR DISPLACEMENT DETERMINATION, 
DATA TRANSFER AND ASSESSMENT. 

10 The present: Invention rela1:es generally 1:o apparatrus for 
and a method of determining a dlsplacemenl: characteristic 
of an object. It also relates to apparatus for and methods 
of transmitting and receiving data* Finally, it relates to 
apparatus for and a method of assessing the approach of an 

15 object to a specified point. 

More particularly, the invention relates to a scoring 
system of the type which can provide a score for the 
proximity of approach of an object such as a missile or 
20 Intruder to a specified location such as a target aircraft 
or other detection location « Although the invention is 
described in detail with reference to a radar system, it 
would also function successfully using, for example, infra- 
red or acoustic techniques. 

25 

In numerous applications a ranging radar system is 
conventionally used to derive measurements of range to one 
or more targets. In some of these applications, the 
accuracy of such range measurements is a principal 

30 requirement. Examples may be found in missile scoring 
radars, collision avoidance and intruder detection. A 
common situation is where a radar sensor is connected to a 
• processor at a ground station via a narrow bandwidth 
telemetry link, and where the sensor must detect the target 

35 with high reliability at near real-time rate. 

An extreme example, which will be used to describe the 
preferred features of the present invention, is that of a 
precision terminal trajectory scoring radar system in which 
40 the differences between range measurements to an object, 
from transceivers separated only by short baselines, are 
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used tro de-termine accurat:ely its angular position, and 
hence provide a proximity score. 

In a conventional radar system, the range of the object is 
5 determined by repeatedly transmitting a radio-frequency 
probe piilse from a transmission means, and receiving this 
at a receiving means. The probe pulse travels both directly 
to the receiving means and also indirectly by way of 
reflection from the object. A detecting means detects all 

10 such signals received by the receiving means, and 
determines the range of the object from the timing delay of 
the reflected pulse relative to the directly received pulse 
(the "direct" pulse). The range is then simply the relative 
delay divided by the speed of light. Such a radar system is 

15 described, for example, in a book by Skolnik, M. I. 
entitled "Introduction to Radar Systems", McGraw-Hill, 
1962. 



In the conventional radar system, the detecting means 
20 samples the signals received by the receiving . means at a 
sampling frequency very significantly less than the 
repetition frequency of the probe pulse. This is 
unavoidable because of various hardware restrictions in the 
system. 

25 

The conventional system suffers from two disadvantages of 
relevance to the present invention. Firstly, because of the 
relatively low sampling frequency, the system cannot cope 
with fast data rates such as might be encountered in 

30 assessing the approach of a hypersonic missile or the like. 
Secondly, the system suffers from what may be termed a 
"dynamic range" problem. The detecting means detects a 
signal which comprises, on the same magnitude scale, both 
the direct pulse and also the reflected pulse. If the 

35 reflected pulse is very much smaller in magnitude than the 
direct pulse, which is frequently the case in systems for 
the detection of small missiles and like small objects. 
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•then the detecting means may have great difficulty in 
accurately detecting both signals, due to their wide 
dynamic range. 

5 A missile scoring radar system is knovm from' International 
Patent Application No- PCT/GB90/00602, which names 
Cambridge Consultants Limited as patent applicant, and 
whose disclosure is incorporated herein by reference. 

10 In this document is disclosed a missile scoring system 
based on an impulse radar technique where the missile is 
sensed as it crosses a sequence of precision range gates 
(or "regions"). From the crossing seqtience of the range 
regions the displacement characteristics of the missile, 

15 such as its trajectory and attitude, can be determined. The 
precision range gating is such that it is possible to 
measure accurately the time at which a feature of the 
approaching missile, such as its nose or tail fins, 
coincides with the range of the range region. Since the 

20 range of the region and the time of coincidence are known, 
they can be used as data from which the position of the 
missile and its trajectory may be calculated. The accuracy 
of the time and range measurements determines that of the 
final trajectory reconstruction. Analogous processes could 

25 be used in other ranging systems. 

In the context of the present invention there are three 
important features of this known system. 

30 Firstly, in order to determine the range of the object, a 
detection timing signal is generated at a delay after the 
time at which the radar probe pulse is transmitted (the 
transmission time), so as to a detect the object when it is 
at a particular range region. That portion of the probe 

35 signal reflected from the object which occurs at the delay 
of the timing signal is then detected, and is used to 
assess the presence or absence of the object at any 
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particular time. A measure of the range of the range 
region^ and hence the object (if present), is then 
determined from the delay of the timing signal. 

5 Whilst this feature has been found to perform well in 
practice, it suffers from a limitation concerning the 
manner in which the range of the range regions is 
calibrated. Calibration is necessary, because in the known 
scoring radar system no absolute determination of range is 

10 made whilst the system is in operation; the delay of the 
timing signal needs to be pre-calibrated since it is not 
measured with reference to any feature of the radar signal 
having a known timing. This is in contrast to the 
conventional radar system described previously above, where 

15 range is determined directly from the time delay between 
the reflected pulse and the direct pulse. 

In the known scoring radar system, the range of the range 
regions is in fact pre-calibrated manually at the 

20 manufacture stage by positioning a reflective sphere at the 
various range regions. This procedure can, however, be both 
cumbersome and inaccurate- Furtheirmore, pursuant to the 
present invention, it has now been fqund that this 
procedure is disadvantageous, in so far as the calibration 

25 cannot be updated during actual operation of the system. 
It has been found that the initial manual calibration may 
drift appreciably due to temperature, ageing and like 
effects in the components of the system. The effect of such 
drift is to introduce inaccuracies into the range 

30 determination. 

Aside from careful manual pre-calibration, and the use of 
crystal oscillators and high-speed logic to determine the 
range region positions, no further special provision is 
35 made in this known system to achieve a high level of 
accuracy for the range measurements. This has been found to 
be adequate when such a system is used on a full-scale 
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1;arget aircraft; however, for use on small subscale tiargets 
(perhaps, for example, aircraft with a wing-span of as 
little as 1 metre ) , greater accuracy may be needed to 
compensate for the smaller baselines offered by such 
5 targets. With smaller baselines/nose to tail dimensions, 
accuracy degrades, since, amongst other things, it becomes 
more difficult to determine accurately the entry angle of 
the missile. 

10 The present invention therefore seeks to provide a means of 
accurately calibrating such a radar system, using such 
range regions. It also seeks to monitor changes in its 
function, since it is important to ascertain the health of 
the system. 

15 

A second significant feature of the known scoring radar 
system is the high rate at which data needs to be 
transferred from the target aircraft to the receiving 
ground station where the data is processed in detail. In 
20 fact, -the known system requires a data transmission 
telemetry bandwidth which is sufficiently wide to preclude 
its use in many potential applications. 

The third and final significant feature of the known system 
25 is the manner in which the approach of the missile is 
assessed. In this known system, the detecting means 
monitors at a high data rate a region awaiting the 
occurrence of an event (in this case the arrival of a 
missile). In practice, the raw data rate from the detecting 
30 means is high (5-20 Mbit/s) and the final data processing 
is complex. Thus it is necessary to identify the occurrence 
of events with a low bandwidth process to avoid overloading 
the complex processor with spurious events. With this known 
system, this has been achieved as follows. 

35 

As explained already, multiple range regions are employed 
to detect the missile when it flies close to the target- 
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When -the missile Is detrected In a range region^ It: Is 
assumed that: larger energy per unit time would be observed 
than when only background noise is detected. Therefore, the 
engagement is identified by summing the detections at all 
5 of the channels (that is, at all the range regions and for 
all of the plurality of receivers which are used). An 
engagement is scored if the sum is above a given threshold. 

Unfortunately, the above procedure has been found to be 
10 insufficiently robust to practical noise mechanisms. In 
particular, in the event of an external interfering 
transmission being received, then typically all the 
channels would detect high signals simultaneously. Also, in 
the event of clutter being received due to reflections of 
15 the radar transmit signal around the aircraft (multi-path 
reflections), again simultaneous events unconnected with 
the approach of the missile would be observed at multiple 
range regions. 

20 In its various aspects, the present invention seeks, 
amongst other things, to provide enhancements to the 
performance of a range measuring/scoring radar system using 
impulse radar techniques. 

25 In overview, these aspects relate to the following 
enhancements . 

In one aspect of the present invention, preferably a real- 
time auto-calibration system is provided which can 
30 continuously measure the range of the range regions and 
hence take out any drift in the system calibration due to 
temperature, ageing or like effects, . 

Secondly, there are requirements where the telemetry link 
35 from the target aircraft to the ground station is of 
limited bandwidth and the data hence has to be compressed. 
Hence in a second aspect there is preferably provided a 
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dat:a compression t:echn±que which can reduce the data by a 
factor of four while maintaining high fidelity and minimum 
latency • In this context, "latency" refers to the delay 
between the generation of the data and Its transmission. 
5 Minimum latency is particularly important in situations 
such as where the target aircraft might be destroyed and 
hence data transmission could end prematurely. 

Thirdly, there are requirements where it is necessary to 
10 produce a near real-time engagement score. Thus in a third 
aspect there is preferably provided an algorithm which can 
enable the missile detection data to be analysed at real- 
time rate using practical computer hardware. 

15 Hence, in one aspect of the present invention, there is 
provided apparatus for determining a displacement 
characteristic of an object, comprising means for 
transmitting, at a given transmission time, a probe signal 
towards the object; means for receiving the probe signal 

20 returned by the object; means for generating a detection 
timing signal at a delay after the transmission time, 
corresponding to a selected range for the object; first 
detecting means, coupled to the receiving means and 
responsive to the timing signal, for detecting the returned 

25 probe signal occurring at the delay of the timing signal; 
and second detecting means for detecting the relative 
timing of the detection timing signal and of a signal 
having a predetermined timing with respect to the 
transmission time; whereby the object displacement 

30 characteristic can be determined from the relative timing 
of these two signals. 

The object displacement characteristic may, for example, be 
its range from a specified location or its velocity. 

35 

In so far as the apparatus of the present invention 
includes a first detecting means to detect (only) the probe 
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signal occurring al: the delay of "the -timing signal, it: 
shares -the two advantages of the system disclosed in 
PCT/GB90/00602 of being able to cope with rapid movements 
of the object and of being able to cope with return signals 
5 of low strength compared to the strength of the direct 
signal. 

However, unlike the known system, a second detecting means 
is provided for detecting the relative timing of the 

10 detection timing signal and of a signal having a 
predetermined timing with respect to the transmission time, 
from which relative timing information a measure of the 
selected range (and hence the object displacement 
characteristic) can be determined. Hence the apparatus of 

15 the present invention can "auto-calibrate", that is, it can 
without manual calibration determine the selected range 
corresponding to the detection timing signal; the apparatus 
can also do this on a constantly updated basis. This can 
avoid the need for an initial manual calibration and can 

20 reduce (or even eliminate) problems due to drift in the 
calibration on account of temperature, ageing or other like 
effects. 

A further advantage which stems from the above is that, 
25 since calibration drift problems can be reduced or 
eliminated, cheaper components can be used, even though 
such components may be prone to significant temperature or 
ageing effects. Hence the cost of the apparatus can be 
reduced significantly. 

30 

Preferably, the transmitting means is adapted to transmit 
the probe signal at regularly repeated transmission times, 
at a given repetition frequency, and the second detecting 
means is adapted to sample the signal of predetermined 
35 timing at a sampling frequency lower than the repetition 
frequency. By arranging that the sampling frequency is 
lower than the repetition frequency, the apparatus of the 
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present ±nvent:±on need not be limited as to the speed of 
object it can detect by the rate at which the calibration 
process takes place. It has been found that in practice, 
due to hardware restrictions, the calibration process 
5 cannot easily be operated at the same rate at which it 
would typically be desirable to operate the repetition 
frequency. Further, it has been found that it is not even 
desirable to operate the calibration process at such a 
rate, since the calibration would typically not need to be 
10 updated particularly frequently. An update rate of once 
every . few seconds or even minutes might typically be 
acceptable . 

Preferably, the sampling frequency is much less than the 
15 repetition frequency, more preferably less than 100 
mlllionths of the repetition frequency, even more 
preferably less than 20 or even 5 mlllionths. In the 
preferred embodiment it is one millionth or less. Lower 
ratios still would be desirable but are not easily achieved 
20 in practice. By operating at such a low relative frequency, 
the requirements on the apparatus in terms of the storage 
and transmission of data can be very significantly reduced. 

Preferably, the apparatus includes two oscillators for 
25 timing the operation of the transmitting and second 
detecting means so that the repetition frequency is timed 
by one oscillator and the sampling frequency is timed by 
the beat frequency between the two oscillators. This is a 
particularly simple and effective way of putting the 
30 invention into practice. 

Preferably, the transmitting means is adapted to transmit 
the probe signal not only to the object but also via a path 
of fixed length, which does not include the object, through 
35 the ambient medium to the receiving means, and the second 
detecting means is adapted to detect the fixed path length 
signal as the signal having a predetermined timing with 
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respect -to -the -transmission -time* By detecting the timing 
of the detection timing signal relative to that of the 
fixed path length signal (the "direct" signal), the 
calibration afforded by the present invention can be 
5 particularly accurate. The "predetermined timing" of the 
direct signal with respect to the transmission time is 
simply the length of the path of fixed length (that is, 
typically the direct distance between the transmitting and 
receiving antennae or else possibly the distance between 
10 the two antennae via a third point such as a reflector), 
divided by the speed of light. This distance can be 
predetermined to a high degree of accuracy. 

In an alternative preferred embodiment, the probe signal or 
15 a transmit trigger signal might, for example, b 
transmitted along a very high quality cable of known length 
and good temperature and ageing properties to the receiving 
means, as the signal having a predetermined timing with 
respect to the transmission time. In this case the 
20 predetermined timing delay of the signal as received by the 
receiving means would be dependent at least in part on the 
length of the cable. 

Preferably, the apparatus includes splitter means for 
25 splitting the signal received by the receiving means to 
both the first and the second detecting means. In the 
preferred embodiment, this signal would include not only 
the returned probe signal but also the direct signal* This 
is a particularly simple way of putting the invention into 
30 effect. 

Preferably, the apparatus includes means for combining the 
signal having a predetermined timing and the detection 
timing signal, and for passing the combined signal to the 
35 second detecting means. 

Preferably, the apparatus includes means for producing a 
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signal represen-tatrlve of t:he stratus of the apparatus In 
dependence on the output of the second detecting means* 
This feature is of particular benefit, since, unlike the 
output from the first detecting means, the output from the 
5 second detecting means, containing as it does information 
concerning the timing of the transmission of the probe 
signal, can provide useful status information, such as 
(perhaps most importantly) whether the transmitting means 
is actually functioning correctly. In the context of 

10 missile tests, status information is particularly important 
in view of the high costs involved in the tests. Very often 
the tests can result in the destruction of the target 
aircraft. If the apparatus of the present invention were 
not functioning properly during such a test, the test would 

15 usually have to be repeated, at great expense. 

Preferably, if the generating means is adapted to generate 
a plurality of detection timing signals (corresponding to 
a plurality of range regions), the second detecting means 
20 is adapted to detect the timing of each of the detection 
timing signals relative to the signal having a 
predetermined timing. This can reduce the amount of 
calibration data produced. 

25 Preferably, the apparatus includes means for determining a 
measure of the selected range from the relative timing of 
the signal having a predetermined timing and the detection 
timing signal. 

30 The invention extends to a method of determining a 
displacement characteristic of an object, comprising 
transmitting, at a given transmission time, a probe signal 
towards the object; receiving the probe signal returned by 
the object; generating a detection timing signal at a delay 

35 after the transmission time, corresponding to a selected 
range for the object; a first detecting step wherein the 
returned probe signal occurring at the delay of the timing 
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signal ±s de-tected; and a second de-tectlng st:ep wherein -the 
rela-tive triming of the detection timing signal and of a 
signal having a predetermined timing with respect to the 
transmission time is detected; whereby the object 
5 displacement characteristic can be determined from the 
relative timing of these two signals. 

Analogous method steps to the above described apparatus 
features are also provided within the scope of the present 
10 invention. 

In a second aspect, the present invention provides 
apparatus for transmitting data. A form of such apparatus 
is known, in which input means is provided for receiving a 

15 data stream comprising a plurality of data samples, and a 
data value (such as a single data bit) characteristic of 
each data sample is evaluated. The data value represents a 
"compressed" form of the data stream. Compression permits 
transmission of the data at a higher rate or over a 

20 narrower telemetry bandwidth. 

The known apparatus has the disadvantage that the 
compressed data values can lag or lead the actual . data 
stream quite significantly in time, especially in 
25 circumstances where the data samples are altering rapidly 
in value. 

The present invention seeks to solve this problem. 

30 In a second aspect of the present invention, there is 
provided apparatus for transmitting data, comprising input 
means for receiving a plurality of data samples; means for 
evaluating a respective block data value characteristic of 
each of a plurality of data blocks, each data block 

35 comprising a plurality of data samples; and means for 
transmitting the respective data value for each block. 
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By arranging the data samples in blocks and transmitrting, 
for each block, a data value characteristic of that block, 
it has been discovered pursuant to the present invention 
that the problem of time lag or lead encountered in the 
5 known data transmission apparatus can largely be obviated* 
The reason for this is that the apparatus of the present 
invention has a longer time (related to the block length) 
over which it can react to changes in the data (especially 
sharp rises or falls ) « 

10 

Each data block may typically comprise a plurality of time- 
sequential data samples. 

Preferably, the evaluating means is adapted to evaluate 
15 individual data values representative of particular data 
samples in each block, and the transmitting means is 
adapted to transmit said individual values. By evaluating 
individual data values representative of some or all of the 
individual data samples and transmitting these, the 
20 accuracy of the technique can be improved still further. 

More preferably, the evaluating means is adapted to 
evaluate the block and individual data values respectively 
as an exponent which is constant for each block and a 

25 plurality of mantissas which are variable within each 
block- This is a particularly effective way of putting the 
invention into practice. It can allow for significant 
variation in overall magnitude between individual blocks, 
whilst also allowing a fair degree of variation actually 

30 within these blocks. 

Preferably, the evaluating means is adapted to evaluate the 
individual data values as taking as possible values the 
value zero as well as at least one other value. This 
35 arrangement can allow data transfer of greater accuracy, 
especially in situations where the stream of data samples 
remains at a largely constant level over a substantial 
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period of t:±me. 

Preferably, the evaluating means is adapted to evaluate the 
individual data values as taking one of an odd number of 
5 possible values. This is especially preferred if one of the 
possible values is zero, so that the compressed 
representations of the data samples can be symmetrical 
about a central zero value. 

10 Preferably, the transmitting means is adapted to transmit 
the individual data values as n bit words, each word 
representing m individual data values, the number of 
possible values for the individual data values taken to the 
power m being no less than 25 % lower than, preferably no 

15 less than 10 % lower than, and no greater than, two taken 
to the nearest power of n. This can ensure that the data 
transfer takes place at a fast and efficient rate, by 
arranging that each transmitted word is as fully packed 
with data as possible. 

20 

In the preferred embodiment, m is 3 , and said number of 
possible values is either 3 or 5, In this event, n would be 
either 5 or 7, so that two to the power n would be either 
32 or 128; the number of possible values taken to the power 
25 m would be either 27 or 125. 

Preferably, the evaluating means is adapted to evaluate the 
individual data values in dependence on the difference 
between data samples. By this feature, the effect of any 

30 d.c- offset in the data stream can be removed. If this 
feature is provided, the transmitting means may be adapted 
to transmit at least one further data value characteristic 
of the absolute value of a given data sample or samples. 
This latter feature can reduce the problems caused by fade 

35 occurring in the data transmission, by providing an 
absolute value from which the receiving apparatus can 
recover the data stream. 
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The inven-tion extends to apparatus for receiving data^ 
comprising means for receiving the data in the form of a 
plurality of block data values, each block data value being 
characteristic of a respective data block, the block 
5 comprising a plurality of data samples; means for 
evaluating the data samples In dependence on the block data 
values; and means for outputtlng the data samples. 

Preferably, the receiving means is adapted to receive 
10 individual data values representative of particular data 
samples in each block, and the evaluating means is adapted 
to evaluate the data samples in further dependence on the 
individual data values . 

15 Preferably, the evaluating means is adapted to evaluate the 
data samples in dependence on block and individual data 
values which are respectively an exponent which is constant 
for each block and a plurality of mantissas which are 
variable within each block. 

20 

Preferably, the evaluating means is adapted to evaluate the 
data samples in dependence on individual data values which 
take as possible values the value zero as well as at least 
one other value. 

25 

Preferably, the evaluating means is adapted to evaluate the 
data samples in dependence on individual data values which 
take one of an odd number of possible values. 

30 Preferably, the receiving means is adapted to receive data 
as n bit words, each word representing m individual data 
values, the number of possible values for the individual 
values taken to the power m being no less than 25 % lower 
than, preferably no less than 10 % lower than, and no 

35 greater than, two taken to the nearest power of n, and is 
further adapted to convert each received word into the m 
individual data values. 
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Preferably, the evaluating means is adapted to evaluate the 
data samples on the basis that the individual data values 
have been evaluated in dependence on the difference between 
data samples. 

5 

Preferably, the receiving means is adapted to receive at 
least one further data value characteristic of the absolute 
value of a given data sample or samples.' 

10 The invention also extends to data transfer apparatus 
comprising data transmitting and data receiving apparatus 
as aforesaid. 

The invention further extends to a method of transmitting 
15 data, comprising receiving a plurality of data samples; 
evaluating a respective block data value characteristic of 
each of a plurality of data blocks, each data block 
comprising a plurality of data samples; and transffldtting 
the respective data value for each block. 

20 

The invention further extends to a method of receiving 
data, comprising receiving the data in the form of a 
plurality of block data values, each block data value being 
characteristic of a respective data block, the block 
25 comprising a plurality of data samples; evaluating the data 
samples in dependence on the block data values; and 
outputting the data samples. 

Method steps analogous to the above described apparatus 
30 features are also provided within the scope of the present 
invention. 

In a third aspect of the present invention, there is 
provided apparatus for assessing the approach of an object 
35 to a specified location, comprising means for detecting the 
likely presence of the object in a plurality of range 
regions at specified ranges from said location, and for 
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producing a detrectrion signal for each range region In which 
t:he likely presence of ^he object: is detected; means for 
determining (preferably Independently) a measure of the 
velocity of the object; and means for assessing the 
5 approach of the object in dependence on whether the 
velocity measure is consistent with the evolution pattern 
of the detection signals. 

By determining a measure of the velocity of the object and 
10 assessing the approach of the object in dependence on 
whether the velocity measure is consistent with the 
evolution pattern of the detection signals, the apparatus 
of the present invention can assess the approach of an 
object with a lower false alarm rate than is possible with 
15 the system known from International Patent Application No. 
PCT/GB9O/00602. The velocity measure can effectively 
provide a way of verifying the presence of the object, even 
in the presence of a substantial amount of noise. 

20 In fact, this aspect of the invention is likely to be 
sufficiently reliable in practice that it can provide a 
fully automatic assessment, without the need for manual 
intervention (as is required with the known system). 

25 Another advantage of the apparatus of the present invention 
is that, as described in detail later, it can operate 
virtually in real-time. This may be important, say, if 
multiple missile firings are taking place, so that the 
results of one firing can be assessed before another 

30 missile is fired. 

The assessment of the approach of the object may, for 
example, and most Importantly, include an assessment of 
whether the detection signals actually do represent the 
35 approach of the object, or whether they represent some 
noise phenomenon or the like. The assessment may also 
include a determination of the closest proximity of the 
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object to the location. 

Means may be provided for producing an assessment signal in 
dependence on the output of the assessing means, 

- 5 

Preferably, the determining means is adapted to determine 
the velocity measure by sensing the velocity of the object 
in a plurality of velocity bands. This can provide an 
apparatus of great versatility, in that it can cope simply 
10 and efficiently with a wide range of velocities. 

The detecting means may be adapted to produce detection 
signals for each velocity band, and means may be provided 
for adjusting the evolution patterns for the velocity bands 

15 so that they are substantially independent of velocity. 
This feature can afford a particularly efficient way of 
assessing the approach of the object, since the evolution 
patterns can then be matched against a single acceptable 
pattern or set of patterns, substantially independent of 

20 velocity* The adjusting means may, for example, operate by 
effecting some form of "downsampling" on the detection 
signals produced for the different velocity bands* 

Preferably, the determining means is such that the velocity 
25 bands overlap one another. By this feature, even if the 
perceived velocity of the object varies somewhat (which is 
likely as it moves past the specified location), it can 
nevertheless be assured that the velocity variation of the 
. object can be accommodated within one such band. 

30 

Preferably, the detecting means is adapted to detect the 
likely presence of the object separately at each velocity 
band, and the assessing means is adapted to ass&ss the 
approach of the object in dependence separately on the 
35 respective evolution pattern of the detection signals at 
each velocity band. By producing such an approach 
assessment, the apparatus can produce more reliable 
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results, and can even assess the approach of multiple 
objects approaching at different velocities. 

Preferably, the assessing means is adapted to assess the 
5 approach of the object in dependence on whether the 
evolution pattern falls between two extremes of 
allowability. This is a particularly simple way of putting 
the invention into effect, and can take account of expected 
relatively minor variations in the evolution pattern from 
10 the theoretical "norm". 

Preferably, the assessing means is adapted to assess the 
approach of the object in dependence on whether the 
evolution pattern approximates to a hyperbolic function. In 
15 so far as a hyperbolic function is a physically realistic 
representation of the evolution pattern (that is, the 
theoretical "norm" ) to be expected in the context, of a 
missile scoring radar system, by this feature can be 
afforded an accurate assessment of the object approach. 

20 

Preferably, the detecting means comprises a plurality of 
individual detecting means, preferably at separate 
locations, and the assessing means is adapted to assess the 
approach of the object in dependence on the respective 
25 evolution pattern of the detection signals for each 
individual detecting means. This feature can ensure an even 
more accurate assessment of the object approach. 

Preferably, the apparatus includes means for generating an 
30 estimate of background noise, and the detecting means is 
adapted to produce the detection signal according to 
whether the detecting means detects the likely presence of 
the object above a noise threshold which is dependent on 
the noise estimate. If so, the generating means may be 
35 adapted to generate a noise estimate which is time 
dependent, so that account can be taken of any changes 
which might occur in the background noise levels. 



wo 94/24580 



PCT/GB94/00738 



-20 - 

The noise tlizreshold may be based on a delayed version of 
the noise estiraatre. Hence the approach of the object need 
not incorrectly influence the noise estimate. For the same 
reason^ the noise threshold may be based on a version of 
5 the noise estimate averaged over a given duration. 

Preferably^ the detecting means is adapted to produce the 
detection signal according to whether the detecting means 
detects the likely presence* of the object above both of two 

10 noise thresholds, one of which is based on the sum of the 
noise estimate and a first constant, the other of which is 
based on a second constant times the noise estimate. This 
has been found to be a particularly Satisfactory way of 
putting the invention into practice, since it takes into 

15 account the different types of noise likely to be 
encountered . 

The invention extends to a method of assessing the approach 
of an object to a specified location, comprising detecting 

20 the likely presence of the object in a plurality of range 
regions at specified ranges from said location, and 
producing a detection signal for each range region in which 
the likely presence of the object is detected; determining 
a measure of the velocity of the object; and assessing the 

25 approach of the object in dependence on whether the 
velocity measure is consistent with the evolution pattern 
of the detection signals. 

Method steps analogous to the above described apparatus 
30 features are ialso provided within the scope of the present 
invention. 

Other preferred features of the invention are as follows^ 

35 In one aspect, the invention consists in a Vernier time 
delay calibration method in which two stable oscillators 
are set at frequencies which differ by a small fraction. 
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One is used 1:o generate a repetitive waveform with at least 
one feature at known delay and other features at unknown 
values of delay (such as in a radar system), whilst the 
other, at a slightly lower (or possibly higher) frequency, 
5 is used to generate a repetitive event whoWe delay with 
respect to the known events of the first clock increases 
regularly from cycle to cycle of the first clock, and by 
coincidence can be used to ' measure the delay of unknown 
events due to the first clock. 

10 

Preferably, the second stable oscillator is derived from 
the first by means of synthesis by a phase-locked loop. 
Alternatively, the second stable oscillator may be 
independent of the first. Again, the second stable 

15 oscillator may be derived from the first by accurate 
control of the beat frequency between the two. Yet again, 
the second stable oscillator is derived from the first by 
accurate control of the interval between coincidences 
between events generated by the second clock and the known 

20 events generated by the first clock. 

Preferably, the "unknown" events are pulses which are 
simultaneously used to close momentarily sampling switches 
which themselves create "range gates" in a radar system. 

25 

In another aspect, the invention consists in a data 
compression scheme where the data is compressed using a 
block floating point scheme where samples are represented 
by a limited length exponent and a limited length mantissa 
30 such that each code represents either the difference 
between successive samples or the absolute value of a 
sample . 

Preferably, the mantissa is a 5 level code with three such 
35 codes packed into a 7 bit word.' 

In another aspect, the invention consists in a signal 
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processor where -the out:pu1: of -the radar Is processed al: a 
discrete number of ranges in a set: of frequency bands 
corresponding to different velocities of the target to 
achieve detection of targets which have a very wide range 
5 of target velocities. 

Preferably, the velocity bins are formed using the Fast 
Fourier Transform and combining the squares of the output 
of the transform in overlapping frequency bands. 
10 Alternatively, the velocity bins could be formed using 
"wavelet transforms", 

Preferably, detection decisions are m^de on each channel 
whenever the signal exceeds an absolute threshold or some 
15 multiple of the noise in each velocity bin. The detection 
threshold may be based on a delayed version of the multiple 
of the noise. 

Preferably, the output of the different-. velocity bins is 
20 down-sampled at different rates to make - the pattern of 
detections velocity independent. 

Preferably, detection decisions are made only as a result 
of a sequence of consistent events across multiple 
25 channels. Consistency here means that for the velocity 
under consideration the sequence of ranges measured is 
consistent with the current velocity. 

Preferably, the final detection decision is made according 
30 to there being consistent events across multiple receivers. 
Here "consistent events across multiple receivers" means 
the same velocity bin at the different receivers at a 
spacing consistent with the spatial separation of the 
receivers . 

35 

The invention extends to a missile scoring system in which 
range calibration is provided by a Vernier time delay 
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calibration method as aforesaid. 

The Invention also extends to a missile scoring system in 
which data compression is achieved according to the data 
5 compression scheme as aforesaid. 

Again, the invention extends to a missile scoring system 
using velocity selective processing according to the signal 
processor as aforesaid. 

10 ) 

Further, the invention extends to a ranging impulse radar 
where calibration is achieved using the venier time delay 
calibration method as aforesaid, where data compression is 
achieved according to the data compression scheme as 

15 aforesaid, or where velocity selective processing is used 
according to the signal processor as aforesaid. 

The various features of the various aspects of the present 
invention may be combined with one another In any 
20 appropriate fashion. 

Preferred features of the present invention will now be 
described, by way of example, with reference to the 
accompanying drawings, in which: 

25 

Figure 1 is a block diagram of a known missile scoring 
radar system; 

Figures 2(a), 2(b), 2(c) and 2(d) (prior art) show, 
respectively, the waveform of a direct wave transmitted by 
30 a first transmitter, the waveform of a direct wave 
transmitted by a second transmitter, and the two switching 
pulses corresponding to the aforementioned direct waves; 

Figure 3 is a block diagram of the scoring radar 
system of the present invention; 
35 Figure 4 is a block diagram Illustrating a data 

compression apparatus according to the present invention; 
Figure 5 is a plot of frequency range against velocity 
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bin, -to illus-tratre the overlap of the Doppler bands in the 
present invention; 

Figure 6 illustrates the best and worst cases of how 
a missile radar signature fits within a time window 
5 employed in the present invention; 

Figure 7 is a plot of radial velocity against scalar 
miss distance for each of six range gates utilised with the 
present invention; 

Figure 8 is a plot of bistatic range against time; 
10 Figure 9 illustrates two acceptable target evolution 

patterns ; 

Figure 10 illustrates first and second points of 
engagement used in the present invention; and 

Figure 11 is a block diagram of a signal processor of 
15 the present invention. 

Different aspects of the present invention will now be 
described in the succeeding numbered sections. The 
description relates specifically to a missile scoring radar 
20 system, but it will be understood that the principles 
adduced can be applied more generally, 

1. aUTO-CAI/IBRflTION 

25 Referring first to Figure 1, a known scoring radar system 
comprises essentially two transmitting means Txl and Tx2, 
two transmitting antenna means 3 and 4, two receiving 
antenna means 5 and 6, two receiving means Rx7 and Rx8 and 
their associated pulse generators PG28 and PG30, clock 

30 means CK9, divider means DVIO, two delay means DLll and 
DL12, a crystal oscillator 13 for controlling the clock 
means CK9 at a constant rate F, and two output means 14 and 
15 for out putting range gate output. 

35 As used herein, and unless the context otherwise demands, 
the terms "transmitter" and "receiver" connote the 
combinations respectively . of the transmitting and 



wo 94/24580 



PCT/GB94/00738 



-25 - 

transml1:t:lng antenna means, and the receiving and receiving 
antenna means. 

The transmitting means Txl, receiving means Rx7 and output 
5 means 14, (together with their .associated circuitry), on 
the one hand, and, on the other hand, the transmitting 
means Tx2, receiving means Rx8 and output means 15 
(together again with their associated circuitry) operate as 
separate first and second sub-systems. Henceforth, 

10 properties relating to the first and second sub-systems are 
denoted by the subscripts "1" and "2" respectively. The 
first and second sub-systems are mounted respectively on 
the upper surface and the lower surface of the target 
aircraft (not shown), so that coverage of the entire range 

15 space around the aircraft is maximised. The relative 
operational timing of the two sub-systems is determined by 
the clock CK9, divider means DVIO and delay means DLll and 
DL12. The transmitting means Txl and Tx2 are adapted to 
provide short radio- frequency pulses when triggered by the 

20 output of the divider means DVIO in the form of transmit 
trigger pulses. The divider means DVIO serves amongst 
other things to reduce the frequency of the oscillator 13 
output . 

25 Typical radio- frequency pulses are shown in Figures 2(a) 
and 2(b) for the first and second sub-systems respectively. 
Such pulses are transmitted from the transmitting means Txl 
and Tx2 via the transmitting antenna means 3 and 4, the 
ambient medium, and thence to the receiving antenna means 

30 5 and 6 and the two receiving means Rx7 and Rx8. A pulse 
which is transmitted directly from the transmitter to the 
receiver antenna means is termed a "direct" wave (or 
"direct" signal), whilst the transmitted pulse as reflected 
off a missile is termed a "return signal". 

35 

The delay means DLll and DL12 are arranged to produce 
detection timing signals in the form of switching pulses 
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respectively at nominal delay times and T2 after 

transmission of the relevant pulse, as shown in Figures 
2(c) and 2(d). The delay means are arranged to trigger 
their respective receiving means Rx7 and Rx8 to act as 
5 switches to couple that portion of the return signal 
occurring when the switching pulses exceed thresholds 
designated (e) and (f) in Figures 2(c) and 2(d) 
respectively into the respective output means 14 and 15. 
In the known scoring radar system, the duration of the 
10 switching Is approximately 200 picoseconds. 

The delays and Xj at which the return signals are output 
to the output means 14 and 15 define two* range gates (range 
regions) at fixed ranges from the target aircraft, the 
15 range being determined from the delay and from the speed of 
light ( "C" ) , The range thickness of the range gates is 
determined from the duration of the switching pulses. The 
basis of the known scoring radar is to detect the presence 
of the incoming missile in the particular range gate. 

20 

In the known scoring radar system, in reality a large 
number of range gates are provided by employing for each 
sub- system a plurality of delay means having different 
delays (in other words, one receiver and transmitter per 

25 several delay means ) , and by employing more than two sub- 
systems, each having its own transmitter and receiver. In 
principle, if the system is to locate the missile precisely 
in 3-dimensional space, a minimum of four sub-systems would 
be required. In practice, eight sub-systems are used, four 

30 for the lower half of the target aircraft, four for the 
upper half, so that a detection can be made regardless of 
whether the missile approaches from below or above the 
aircraft. 

35 In the known scoring radar system, cabling (for instance 
between the transmitting means Txl,Tx2 and the transmitting 
antenna means 3,4), the electronic transmitting elements 
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( for instance in the transmitting means ) , the receiving 
amplifiers in the receiving means and printed wiring tracks 
all introduce nominal transmit delays into the two sub- 
systems, denoted respectively St^ and etj. The delays are 
5 due, for example, to charging delays owing to the 
capacitance of the cabling and the electronic transmitting 
elements. The delays are subject to small unknown timing 
errors Ct^ and etj, attributable to temperature, ageing and 
geometric variations and like effects. The pairs of 
10 transmitting and receiving antenna means 3, 4, 5, 6 are 
separated by accurately-known distances d^ and dj, so that 
the signals arriving at the receiving means Rx7, Rx8, 
directly from the transmitter means Txl, Tx2 do so at 
delays t^ and tj given by 

15 

tj « 6ti + €ti + dj/C and ~ et, + etj + dj/C 

Likewise, the printed wiring tracks and electronic devices 
in the receiving means Rx7, Rx8 introduce receive delays dr^ 
20 + and drj + Erj into the timing of the detection timing 

signal. Thus the actual delays of the detection timing 
signals (switching pulses) are 

+ 6rj + €rj, and Xj + Srj + Erj 

25 ■ 

The range gates referred to previously are in the form of 
ellipsoids with foci at the respective transmitting and 
receiving antenna 'means 3, A, 5, 6 such that the sum of the 
distances from the transmitting antenna means to any point 
30 on the surface of the missile and back to the receiving 
antenna means is equal to R, given, for the first and 
second sub -systems respectively by: 

Ri = C * (Xj + 6r, + er^ - ( dt^ + ^tj)), and Rj* C * (X2+ ^^r^ 
35 + erj - (6t2 + etj) ), 



where C is the speed of light, as stated above. 
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In -the knovm system, the values of and are measured by 
a one-off manual calibration^ leading to estimates of (er^ - 
6ti ) and ( 6T2 - 6t2 ) . The uncertainty in timing arising from 
thermal and cable variations and the like is then er^^-et^ 
5 and erj-Gt^, and the differential timing uncertainty between 
the two range gates of the first and second sub-systems is 
(6ri - €1, - €ti + etj), 

A preferred .embodiment of the present invention is now 
10 described with reference to Figure 3. In Figure 3, the 
same components as are shown in Figure 1 are represented by 
the same reference numerals. Only those components 
additional to those described with reference to Figure 1 
are now described in detail. 

15 

The improvement which is the subject of a preferred aspect 
of this invention is to introduce a third and fourth 
receiving means Rxl9 and Rx20, splitting means 16 and 17, 
a further clock means 24, a further crystal oscillator 18, 

20 further dividing means DV25 and DV23, further pulse 
generation means PG21 and PG22, attenuation means 33 and 
32, adding means 29 and 31 and further output means 26 and 
27. The further clock means 24 is controlled at a rate Fj, 
where Fj « F^-f . F^ is the frequency of the clock means CK9, 

25 and f is also a constant frequency which may be a small 
submultiple of F, Usually f will be much smaller than F^. 
In the preferred embodiment f is IHz, whilst F^ is of the 
order of 2MHz. 

30 As described above in relation to the r^celvtzig means Rx7 
and Rx8, the further receiving means Rxl9 and Rx20 are 
arranged to receive switching pulses (detection timing 
signals) from the pulse generation means PG28 and PG30, in 
this event via the attenua-tlng means 33 and 32 and adding 

35 means 29 and 31. The receiving means are arranged to be 
switched by the further pulse generation means PG21 and 
PG22 so as to sample at a data rate determined by the 
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fur-ther clock means 24. 

In addition to closing the receiving switches in the 
receiving means Rx7 and Rx8, the pulses (a) and (b) shown 
5 in Figure 2 are in parallel used to make precision timing 
measurements • 

In one preferred aspect, the effectiveness of the 
improvement provided by the present invention resides in 
10 the fact that no significant uncertainties in delay occur 
between the times when each of the signals (c) and (d) in 
Figure 2 exceed the respective thresholds and the 
respective times when the switches in the receiving means 
Rx7 and Rx8 are closed. Any time differences are the 
15 result of fixed and well-known propagation delays in short 
cables or printed wiring within thie electronics unit. Thus 
a measurement of the instant when the respective signals 
(c) and (d) reach a peak or the centre time of the interval 
over which they exceed a threshold can be taken as an 
20 accurate measurement of the instant when the switches in 
the receiving means Rx7 and Rx8 are closed. Such a 
measurement is made in the third and fourth receiving means 
Rxl9 and Rx20. The measurement, if made at intervals of 
teal, provides direct measurements of (6r^ + er^) and (Srj 
25 + Srj), rather than merely and 6r2, increasing the 

accuracy of the range gate calibration, provided that Sr^ 
and 6r2 are stable over the time teal. The same mechanism 
can be used to measure the arrival time of the direct 
transmitted wave at the receiver, giving the values of (6ti 
30 + et^) and (et^ + et^). 

The preferred embodiment of scoring radar system according 
to the present invention is now explained in more detail 
with reference to Figure 3. The system can be understood 
35 with reference to the top half of Figure 3 because the 
bottom half is merely another channel disclosing a second 
transmission to be calibrated for a multistatic system. 
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The transmi-tter Txl pulses regularly in response to -trigger 
pulses from the clock means CK9 and dividing means DVIO, 
The -trigger pulse is delayed by the delay means DLll and 
triggers the pulse generator PG28 to produce a switching 
5 pulse (detection timing signal) at the delay defining the 
range gate position. A small fraction (perhaps 1/lOOOth) 
of the switching pulse defining the delay of the switching 
range gate is fed into the adding means 29 which combines 
a small fraction of the signal from the splitting means 16 

10 and passes this to the third receiving means Rxl9 which is 
clocked at a small offset frequency to the master clock Fj. 
Thus the third receiving means will output a signal to the 
calibration output means 26 which contains the direct 
signal transmitted directly between the transmitting and 

15 receiving antenna means 3 and 5 and the detection timing 
signal corresponding to the given range gate- Hence the 
calibration of the range gate can be derived from the time 
delay between the direct signal and the detection timing 
signal . 

20 

Because the swltch&s in the third and fourth receiving 
means Rxl9 and Rx20 are repeatedly closed at a constant 
rate T^, slightly less than F^, the trigger time of these 
receiving means will slip with respect to the trigger time 

25 of the fixed delay receiving means Rx7 and Rx8. The 
switches in the former receiving means can be regarded as 
samplers, with their outputs connected to the output means 
26 and 27. As this delay (that is, timing slip) increases, 
at one instant the sampling pulse of the receiving means 

30 Rxl9 will coincide with the direct wave (a) (see Figure 2); 
then some time later ("tcgate") it will co-incide with the 
detection timing signal (c), and then at a time ("teal") 
after its coincidence with the direct wave it will co- 
incide again with the direct wave. If the output means 26 

35 and 27 take the form of low-pass filters, the signal 
observed there will be the sum of (a) and (c) in Figure 2. 
Provided that the difference frequency f is constant over 
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the interval teal, the range gate delays (Ri/C) and (Rj/C) 
may be measured accurately, provided also that d^ and d, are 
known equally accurately. Under these circumstances, 

5 Ri/C - dj/C « tcgate/tcal 

1/f 

The analogous expression applies to the lower sub-system 
10 .("2"). 

Note that in this expression the values 6t^, Stj, 6r^, 6r2, 
etj, etj, €ri and Erj do not appear; for. a given tcgate the 
value of Rj depends only on the accuracy of d^, which can be 
15 measured with a tolerance of millimetres, and f , which is 
crystal-controlled. 

This expression is precise provided that both F2 and F^^ (and 
hence f ) are constant over the period teal , and yields 
20 precise, absolute calibration for the two range gates 
described • 



A differential calibration can also be obtained by 
calculating the difference between the range values and 
25 R2 for physically close sub-systems respectively on the 
upper and lower surfaces of the target aircraft. This is 
given by 



(Ri - R2)/C - (di - d2)/C = 6t,^ - 6t„ 

30 

1/f * teal 

In this equation, 6t^^ is the difference between the arrival 
times of the direct waves from the upper and lower sub- 
35 systems (as detected by the third and fourth receiving 
means Rxl9 and Rx20), and 6tg^ is the time difference 
between the detection timing pulses for the upper and lower 
sub-systems. For this differential calibration to be 
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accurate, t:lie strabilitry of the difference frequency f is 
only important: over the shorter time et^t minus 6t„- This 
differential calibration of the range gates may be 
particularly useful in circumstances where the transceivers 
5 on the upper and lower sub-systems are quite close, and 
hence the scoring radar system is quite sensitive to the 
relative, positions of the range gates for the two sub- 
systems. In such a circumstance the standard calibration 
technique might be employed to calibrate the range for a 
10 range gate on one of- the sub-systems, and the differential 
calibration technique might be employed to calibrate the 
range for the corresponding range gate on the other of the 
sub-systems . 

15 The calibration technique described above may be termed a 
Vernier time delay calibration, and in the example is 
designed to provide delay calibration accurate to a few 
tens of picoseconds* 

20 In practice,- the calibration output from the output means 
26 and 27 (Figure 3) is sent to the ground station as a 
separate output from the range gate output from the output 
means 14 and 15. The calibration of the system is updated 
roughly once every eight seconds, which is the rate at 

25 which information from all eight receivers is sampled. The 
system is calibrated according to either the last 
calibration signal or an average of the last few 
calibration signals. 

30 Alternatively, the calibration output may be processed on 
the target aircraft. However, since the calibration output 
may contain useful information about the status of the 
system, processing is preferably carried out in the ground 
station. 

35 

In practice,* the calibration output from the system is low- 
pass filtered and then digitised, say, at a rate of 5kHz. 
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Sampling could be at: a lower rate stdLll If desired. 

As described above, the calibration Is carried out by 
circuits forming an Integral part of the scoring radar 
5 system. These allow continuous calibration "to be carried 
out with a repeat cycle of 1/f. However, the calibration 
could alternatively be carried out using separate 
calibration equipment, whenever convenient. In this latter 
alternative, compensation could not be made for thermal and 
10 like effects between calibrations. 

In the preferred embodiment. It is Important to ensure that 
there is no cross -coupling between the receiving means Rxl9 
and Rx20 (and their associated circuitry) and the receiving 

15 means Rx7 and Rx8 (and their associated circuitry). Such 
cross -coupling might otherwise interfere with the signals 
detected by the two sets of receiving means. Cross- 
coupling can be reduced by ensuring that the number of 
possible conducted or radiated communication paths is kept 

20 to a minimum. 

In the most preferred embodiment, there are in total 
actually 8 detection receiving means (such as Rx7 and RxB) 
provided on the target aircraft. For each of these 
25 receiving means, there are, further, six separate delay 
means ( corresponding to six range gates ) , so that each 
receiving means receives signals from the six separate 
delay means. Hence, there are 48 channels of range gate 
data altogether. 

30 

For each detection receiving means such as Rx7 and Rx8, 
there is one calibration receiving means such as Rxl9 and 
Rx20, with the detection timing signals of each of the six 
delay means being passed to the single calibration 
35 receiving means. Hence there are in total eight channels 
of calibration output. 
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The various channels of data are -transmi-ttied via -the 
•teleme-try link to the ground station in bursts of data. In 
order to avoid unnecessary signal drift, the eight 
calibration outputs are transmitted in bursts in pairs, one 
5 member of each pair relating to a receiving means on the 
lower part of the target aircraft, the other member 
relating to the most closely coupled receiving means on the 
upper part of the target aircraft, and so on. Each burst 
of calibration output contains a waveform of duration 
10 slightly greater than 1/f, so that the waveform contains 
two direct transmit signals. This makes processing easier. 

2. DATA COMPRESSION 

15 The preferred embodiment of missile scoring radar system 
produces a large amount of data which has to be transmitted 
by radio telemetry from the target aircraft to the 
receiving ground station. There are^ simultaneously, data 
from -48 range gate channels (8 receivers each having 6 

20 range gates). The data rate is high, in order to cope with 
fast-moving missiles. The amount of data to be transmitted 
hence requires a large telemetry bandwidth. 

In the preferred embodiment of the second aspect of the 
25 invention, concerned with data compression, a scheme is 
presented for compressing the data during transmission from 
the target aircraft in order to reduce the telemetry 
bandwidth. The data is then decompressed at the ground 
station. Whilst this second aspect is described herein in 
30 the context of a scoring radar system, it is not limited to 
such a system. 

In overview, the data compression procedure is as follows. 

35 1) The analogue signal at each instant is digitised at 
high dynamic range ( 8 bit or greater in the preferred 
embodiment). The sequence of data samples thus digitised 
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is referred to as x^. 

2 ) The digitrised dat:a samples are buffered up into blocks 
of multiple samples. 

5 

3) For each channel and each block, the maximum absolute 
difference between two successive samples is evaluated 
(that is, max Ix^-x^.^^l is evaluated) for all but the first 
sample of the block (see Step 8 below). 

* 10 } 

4) One exponent exp per block (a block data value) is 
defined which when 2 is raised to the power of the exponent 
and then multiplied by (2,1,0,-1 or -2) will exceed the 
maximum difference (max jx^ - x^^^ \ ) evaluated in Step 3 

15 above- 

5) Each data sample is encoded to form an individual data 
value by finding the best value of r from the set (2,1,0,-1 
and -2) which when multiplied by 2"^ minimises the 

20 difference between the current sample and the last 
reconstructed value (see Step 6 below). The sequence of 
encoded individual data values is known as r^. 

6) The new "last reconstructed value" is evaluated as the 
25 previous last reconstructed value plus r times 2"'', and is 

known as y^. 

7 ) Each set of 3 consecutive values of r is combined into 
a respective 7 bit word equal to the sum of (r^.j^^) + 

30 5*(r^.i+2) + 25*(ri+2). 

8) For each new block the first difference used in 
defining the exponent is x^ (that is, from the new block) 
minus y^.j (that is, from the previous block). This ensures 

35 that any residual error from the previous block does not 
accumulate from block to block. 
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The scheme outillned above is herein referred tro as the 
"differential block floating point scheme" • Significant 
features of the differential block floating point scheme 
are now described in greater detail. 

5 

Firstly, the scheme is herein termed a differential scheme, 
hy which, is meant that each compressed, encoded data value 
is evaluated in dependence upon the difference between two 
data samples to be compressed, rather than representing the 
10 absolute value of a data sample to be compressed. This 
allows compensation to be made automatically for any d.c. 
offsets in the data (such as frequently occur in data from 
radar systems ) . 

15 Secondly, the scheme- operates in blocks, with a data ^value 
characteristic of a whole block being evaluated, as well as 
individual data values characteristic of the individual 
data samples. This affords a scheme which reacts quickly 
to sharp variations in the magnitude of the underlying data 

20 stream, and hence follows . the data stream without 
significant time lag or lead. 

In the preferred embodiment, there are either 6 or 9 data 
samples per block. At a data rate of 31kHz a block length 

25 of 6 samples has been found to operate satisfactorily in 
practice without undue distortion. In different 

circumstances, other block lengths might be appropriate, 
although in the preferred embodiment they would need to be 
a multiple of three samples. With short block lengths, the 

30 data samples can be followed more accurately, but the data 
rate is higher. 

Thirdly, the scheme is termed a "floating point" scheme, by 
which is meant that the data is represented in mantissa and 
35 exponent form. 

The exponent (exp) is a data value characteristic of the 
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en-tire block. I*t ±s lnt:ended "to be representatilve of -the 
dynamic range of the data in the appropriate block. In the 
preferred embodiment, the exponent is 3 bits long^ so that 
can vary from 1 to 2' (128). 

5 

In the preferred embodiment, the exponent is the nearest 
integral value of log^ (4/3* maximum absolute difference), 
since this has been found to be optimal. However, for 
different applications, logarithms to different bases might 
10 be appropriate. Alternatively, the exponent might be 
chosen from a look-up table. 

The decoding process, which would typically be carried out 
at the ground station, is the inverse of the above encoding 
15 process. 

Another realisation of the data compression scheme is 
exactly as above except that the exponents are defined from 
the biggest absolute value of the digitised samples x^ 
20 rather than the differences Xi-x^.i- This scheme herein is 
referred to as the "absolute block floating point scheme". 

Returning again to the differential block floating point 
scheme, the mantissa (r) is a (compressed) data value, one 

25 for each individual sample. The mantissa is in fact highly 
compressed, taking one of a set of only five possible 
values, -2, -1, 0, +1, or +2. The existence in the set of 
the value "zero" is significant, since the scheme can 
thereby cope easily with long periods when the data stream 

30 does not vary significantly in value. 

Fourthly, the data compression scheme stores the individual 
data values representative of the individual data samples 
in sets of 3, two sets to a block, as a single 7 bit word. 
35 A 7 bit word can take one of 127 possible different values. 
Since each data value has a possible 5 values, a set of 3 
data values has 5^ (equals 125) possible different values. 
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Hence "the division of -the block into sets in this way can 
permit: economical datra transfer. 

Fifthly, an important feature of the data compression 
5 scheme is its ability to cope with fading of the 
transmitted signal. If the signal fades, clearly it may be 
impossible for the ground station to receive any reliable 
data during the period of fading. This presents a 
particular problem when the differential scheme is used, 

10 since in this event recovery of the true absolute level of 
the data stream is impossible without further information - 
To combat this, every, say, 24 blocks the absolute value of 
a data sample is transmitted, along with a synchronisation 
word to allow the data stream to be re-synchronised. Given 

15 the absolute value, ^ this can then be propagated back in 
time to retrieve data samples back to the point at which 
the transmission fade ceased. 

Sixthly, further data compression may be achieved by using 
20 the same exponent for the— output from ^ some or- all of the 
receiving means. 

The preferred embodiment of data compression apparatus 
according to the present invention is explained with 
25 reference to Figure 4. In the example given there are 48 
channels and the exponent is defined with respect to blocks 
of nine samples, but the radar system could of course have 
any number of channels and the exponent could be defined 
with respect to an arbitrary number of samples. 

30 

The preferred embodiment of data compression apparatus 
comprises a first Random Access Memory (RAM) 201, a second 
RAM 202, a third RAM 203, a first subtracter 204, a fir^t 
Programmable Read Only Memory (PROM) 205, a second 
35 subtracter . 20.6, a second PROM 207, and an adder 208. The 
input to the apparatus is a new data input sample which is 
provided to the line marked "Input"; the output is an 
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encoded (compressed) datza value including bo-th a mantissa 
and an exponent (see the line marked "Output (Mantissa and 
Exponent ) " ) . 

5 In overview, the apparatus operates in the following 
fashion. The new input sample is subtracted from the last 
previous input sample, which has been stored in the RAM 
201, by the subtracter 204. The exponent corresponding to 
the result of this subtraction is compared by the PROM 205 

10 with the current exponent and is updated if the new 
difference is larger. The new exponent is written back into 
the RAM 202 . The corresponding sample in the previous block 
is subtracted by the subtracter 206 from the last 
reconstructed value held in the RAM 203. The resultant 

15 difference is then encoded by the PROM 207 which outputs 
the mantissa and exponent as the output of the apparatus. 
The value of r x 2"** is also output from the PROM 207 to the 
adder 208. This adder adds this value to the last 
reconstructed value and hence updates the last 

20 reconstructed value in the RAM 203. 

The following is a detailed description of the various 
operations which are required to compress 9 data samples 
into a single data block. The description proceeds with 
25 reference to three time-consecutive data blocks, n-1, n and 
n+1. The entire sequence of operations is repeated for 
each of the 48 data channels. 

30 For each of the samples 2 to 9 in block n(n,2 to n,9): 

Increment the sample number from 2 to 9 and the channel 
number from 1 to 48 for each of Stages 1 to 3. That is:- 

35 For sample = 2 to 9 
For channel = 1 to 48 
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STAGE 1 

For block n 

5 Read the previous value for this channel from Random Access 
Memory 201 (n, saraple-1) and latch. 

Read the current exponent from RAM 202 ("new" exponent) and 
latch. 

10 

For block (n-1) 

Read the exponent for the previous block from RAM 202 
("old" exponent) and latch. 

15 

Read the "last reconstructed value" for the previous block 
from RAM 203 (n-1, sample-1) and latch. 

STAGE 2 , 

20 

For block (n-1) 

Read the next sample from the previous block from RAM 201 
(n-1, sample) and latch. 

25 

Calculate the difference between the last reconstructed 
value (n-1, sample-1) and the next sample (n-1, sample) in 
subtracter 206. 

30 Code this difference in PROM 207, and pass the coded word 
(r) to a framing circuit (not shown), which arranges the 
coded words in the desired order for transmission. 

Add the transmitted difference (r x 2"^) using the adder 208 
35 to the. last reconstructed value (n-1, sample-1) to give the 
new "last reconstructed value" (n-1, sample). 
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STAGE 3 
For block n 

5 Wrlt:e a new sample (n, sample) Into Random Access Memory 
201. 

Subtract: the previous value (n, sample- 1 ) using subtracter 
204 and update the "new" exponent using the PROM 205. 

10 

Write the "new" exponent to the RAM 202. 
For block f n-1 ) 

15 Write back the "old" exponent (unchanged) to the RAM 202. 

Write the new "last reconstructed value" (n-1, sample) to 
the RAM 203, 

20 Next channel 

Next sample 

Thus, at this time, one of the RAM's (the RAM 201) contains 
25 the 48 channels of the nth block of 9 samples. We have also 
stored the last reconstructed value of block n-1 and the 
exponents for block n-1 and block n based on all 
differences except the difference between the first sample 
of the nth block and the last reconstructed value of block 
30 n-1. The sequence of operation then proceeds as follows. 

For the first sample in the next block (n-i-l^l): 

Increment the channel numbers from 1 to 48 for each of the 
35 following stages. That is:- 

For channel * 1 to 48 
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STAGE 4 
For block n 

5 Read the "new" exponent from the RAM 202 (for block n) and 
latch. 

Pass the "new" exponent to the coding PROM 207 to form part 
of the coding equation. 

10 

Read the last reconstructed value (n-1,9) from the RAM 203 
and latch. 

STAGE 5 

15 

For block n 

Read the first sample of block n (n,l) from the RAM 201 and 
latch. 

20 

Subtract the last reconstructed value (n-1,9) from sample 
(n,l) using the subtracter 206. Increase the size of the 
"new" exponent (exp) if necessary. 

25 Code the difference (r) using the PROM 207, and pass the 
coded word (r) to the framing circuit. 

Add r X 2'*P using the adder 208 to the last reconstructed 
value (n-1, 9) to give the new "last reconstructed value" 
30 (n,l). 

STAGE 6 

For block f nn-l) 

35 

Write the new sample (n+1,1) into the RAM 201. 
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For block n 

Write the "new" exponent updated in Stage 5 into the "old" 
exponent memory (RAM 202). Write zero into the "new" 
5 exponent memory. * ' / . 

Write the new "last reconstructed value" (n^l) also 
calculated in Stage 5 to the RAM 203. 

10 Next channel. > 

It will be appreciated that the data compression apparatus 
described above could be implemented in several alternative 
15 forms. It might, for example/ be implemented on a 
programmable processor, in discrete digital logic, in an 
Application Specific Integrated Circuit, or in field 
programmable gate arrays . 

20 ) 

3. REAL TIME SIGNAL PROCESSING FOR THE SCORING RADAR SYSTEM 

3.1 Overview 

25 The operation of the scoring radar system known from 
International Patent Application No. PCT/GB90/00602 has 
already been described. In the preferred embodiment of the 
present invention, two main enhancements have been made to 
this system to increase the reliability of the scoring 

30 process, and also to enable signal processing to occur in 
real time . 

The first enhancement is that the energy bandwidth used for 
basic detection is rendered velocity dependent. In other 
35 words, a missile proximity score is derived at each channel 
by measuring the energy in a Doppler frequency bandwidth 
corresponding to the missile velocity. This velocity, which 
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is unknown at "the trime o£ the test:, is estimated by 
measuring the Doppler shift. 

The second enhancement is that, because the range gates of 
5 the scoring radar system are of known spacing, given the 
missile velocity, a prediction can be made as to when a 
missile of a velocity measured at one range gate would be 
expected to arrive at the next range gate. 

10 Hence, for a given receiver, the likely presence of a 
missile can be identified based on a positive detection at 
the same velocity (or related pattern of velocities if 
acceleration is significant) over a number of range gates. 

15 Finally, an* engagement is scored if multiple" receivers 
identify the likely presence of a missile at the same (or 
similar) velocity. 

The multi-stage approach, described above . can afford 
20 significant protection against false alarms- ' caused by 
received interference or clutter, because these phenomena 
are unlikely to be correlated in frequency (and hence 
apparent velocity) at different channels. Furthermore, to 
the extent that noise occurs at different channels at the 
25 same frequency, it is unlikely that it will occur in the 
correct sequence to pass the test of being consistent at 
a given receiver. In the unlikely event that this occurs 
then it is unlikely that the interference will present this 
consistently to the other receivers. 

30 

3*2 Background to the Preferred Embodiment of Scoring 
Radar System 

In the preferred embodiment of scoring radar system, eight 
35 receivers and six range gates are provided, which can 
detect missiles with speeds from 200 to 8000 feet per 
second (roughly 60 to 2500 metres per second). The range 
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gatres are set at 40, 60, 80, 100, 110 and 125 feet (roughly 
12, 20, 25, 30, 35 and 40 metres). Each range gate is a 
mere roughly 1 inch (2.5cm) in thickness. Each gate is 
sampled at a data rate of approximately 35 kHz and the 
5 samples are converted to 10-bit digits. In alternative 
embodiments, the range gate spacing might be as great, say, 
as 100 or even 500 feet (30 or 150 metres) or as little as, 
say, 5 or even 1 foot (1.5 or 0.3 metres). 

10 For data compression purposes (for which see Section 2 
above), the data is encoded using a differential block 
floating point scheme. Each sample has a five level 
mantissa and 3 successive samples are packed into 7 bits 
and send a new exponent for every group of 9 (or possibly 

15 6) samples. 

Given the range of target speeds to be detected and the 
setting of the range gates, the fastest time between range 
gate crossings is 2.5 ms, while the slowest is roughly 
20 100ms. The maximum total encounter time is 1.25 seconds. 

3.3 The Scoring Algorithm 

3.3.1 Overview 

25 

The scoring algorithm described in this section is required 
to perform well in the presence of various types of noise. 
This includes white noise, impulsive noise and narrow band 
interference. To perform well, the algorithm should not 
30 fail to detect a genuine engagement, but should have a low 
false alarm rate. 

The algorithm described herein is applicable to detection 
of a moving object (such as a missile). In the case of 
35 missiles, the relative trajectory of the missile relative 
to the target aircraft is approximately linear and hence 
the range-time curve with respect to an individual receiver 
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(assximlng bo-fch approach and recession of "the missile) is 
approximatiely hyperbolic. The present: algorit:hm functions 
-to mat:ch the detection process to the known geometric 
(approximately hyperbolic) evolution pattern of the missile 
5 with respect to the receivers « This involves looking for a 
sequential pattern of range gate crossings representing the 
approach of the missile. The algorithm also allows the 
possibility of looking for a similar pattern as the missile 
recedes. This feature can be used to ascertain physically 
10 possible detection sequences. 

In addition, as described previously, the algorithm takes 
advantage of the fact that the signals 'of interest from a 
missile travelling at a particular speed are confined to a 
15 limited bandwidthi (This follows from a consideration of 
the Doppler effect. ) 

One of the difficulties to be overcome is the wide range of 
possible missile speeds. This is dealt with by searching in 

20 the frequency domain for cumulative activity corresponding 
to a set of assumed missile speeds which ^ cover the full 
range from 200 to 8000 ft/s (roughly 60 to 2500 m/s). For 
each speed in the set, a separate attempt is made to 
determine whether there is a corresponding set of range 

25 gate crossings. 

The search for sequential patterns of range gate crossings 
(so-called "pattern matching") is carried out independently 
for each receiver. This yields a set of scores, one per 
30 receiver, one per asstimed speed, which are related to the 
number of range gates that have been crossed at that speed. 

The task of amalgamating the scores produced by each 
receiver into an overall assessment of the engagement and 
35 hence an overall engagement score involves looking for 
consistent evidence at the same speed of a likely 
engagement from two or more receivers. 
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In out:llne, 1:he s1:eps performed In -the algorithm are:* 

1) Decompress the data, by inverting the five level 
differential block floating point coding described in 

5 Section 2 above. 

2) Buffer up data in buffer blocks of 512 samples per 
channel. Thus a storage capacity of 24k words is 
required at this point. 

10 ) 

3) Implement Fast Fourier Transforms on the data (without 
overlapping). The FFT processing is invoked twice on 
each block of 512 samples, once tb do 2 x 256 point 
EFT'S and once to do 8 x 64 point FFT's. 
(Alternatively, equivalent processing could be 
accomplished using wavelet transforms, where the 
frequency bins are formed by recursive filters 
implemented in the time domain. } 

Square the output from the FFT processing to obtain 
power spectral estimates. 



15 



20 4) 



5} Sum the spectral estimates to estimate the power in 
various velocity bins, that is, those corresponding to 
25 a series of different missile velocities. Downsample 

these estimates by various amounts (that is, sample 
them at different rates) according to velocity. 

6) Form a background noise estimate. The estimate is 
30 delayed before being used so that recede events can be 

detected. 



7 ) Threshold current energy in the velocity bin bandwidth 
against a multiple of the (background) noise plus a 
35 quantisation limit (or an absolute value if 

appropriate ) . 
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8) Thus ob-tain a detectrion decision for each range ga-te, 
for each velocity bin and each receiver, at a rate 
which is scaled approximately according to velocity 
{ see Step 5 above ) . 

9) For each receiver, look for patterns of detections 
across all range gates that match those expected for 
approach or recede events, to derive an overall 
detection decision, and hence an overall score, for 
that receiver . 

10) Derive an overall engagement score based on the 
overall detection decisions for each receiver. 

15 The main features of these steps are described^ in the 
succeeding sub- section. Apparatus to put the algorithm into 
effect is described in Section 4 below. 

3.3.2 The main processing steps of the scoring algorithm 

20 

For convenience, the scoring algorithm is described as a 
number of processing steps, although there is nothing 
fundamental about the manner in which the steps have been 
divided . 

25 

Step A: Divide signals into frequency bands 

Signal to noise ratio improvement and velocity selectivity 
is gained by dividing the signal into a number of frequency 
30 bands ("bins'.*) (in the present embodiment sixteen). 
Subsequent processing steps will search for a missile of a 
particular speed only in the frequency band in which its 
signature is expected to fall. 

35 The incoming signals are transformed twice by Fast Fourier 
Transform, once in blocks of 64 and once in blocks of 256, 
This is done to achieve a good compromise between time 
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resolution and frequency resolutiion. In the lowest velocity 
band, a longer FFT is needed in order to obtain sufficient 
frequency resolution. 256 samples gives 137 Hz resolution, 
which corresponds to a speed resolution of around 40 m/s 
5 (compared with the slowest possible missile speed of 
roughly 60 ro/s). In the higher velocity band, a shorter FFT 
is needed in order to achieve adequate time resolution. 64 
samples occupies a time of 1.8 ms (at 35 kHz sampling 
frequency), which is comparable with the duration of the 
10 missile signature at the fastest speed. 

As mentioned above ^ an equivalent process could use wavelet 
transforms, to obtain power estimates in the various 
frequency bands* 

15 

The data is transformed to the frequency domain and power 
spectral estimates computed by squaring and adding the real 
and imaginary outputs from the FFT. 

20 The data is processed in 16 velocity bins and thus the 
velocity resolution of the system is ±12.5 %. Other 
numbers of velocity bins may be appropriate, maybe as low 
as 4 or as high as 32 or 64. For the efficiency of the FFT 
processor, a number of bins which is a power of two is 

25 preferable. The power in each velocity bin is computed from 
the sums of the FFT spectral estimates. These velocity 
bin estimates are then downsampled by summing together the 
powers from consecutive time samples. The downsampling 
rate is chosen according to velocity bin such that the time 

30 between samples is comparable with the duration of the 
missile signature at that velocity. This ensures that the 
signal to noise ratio is optimised, and also' permits the 
pattern of movement of the missile to be matched against a 
single pattern, regardless of velocity. Also, for a given 

35 velocity bin, the downsampling rate is preferably chosen 
such that the number of detection samples per range gate 
crossing is in the range 1-2. 
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The "table below shows how the full range of velocities can 
be divided into 16 velocity bins. The table provides, for 
each velocity bin, the median velocity, the bounds of the 
Doppler frequency bands covered by the relevant bin, the 
5 engagement time, -the size of the Fast Fourier Transform, 
the number of FFT bins employed, the downsample rate, and 
the overall downsample rate (computed as the downsample 
rate times the FFT size). 
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Median 

Velocity 

(m/s) 



Doppler Engagement FFT 
(Hz) Time (sec) size 



ffTbins Down- 
sample 
, rate 



Overall 

downsample 

rate 





88 


136, 528 


1.25 


256 


2, 5 


6 


1536 




110 


220, 660 


0.692 


256 


3, 6 


5 


1280 




138 


276,828 


0J52 


256 


4. 8 


4 


1024 




172 


344. 1032 


0.443 


256 


4. 9 


3 


768 


10 


215 


430. 1290 


0J54 


256 


5, 11 


2 


512 




269 


538. 1604 


0.283 


256 


5, 13 


2 


512 




336 


672, 2016 


0.226 


256 


6, 16 


2 • 


512 




420 


840. 2520 


0.181 


256 


8. 20 


1 


256 




524 


1048, 3144 


0.145 


64 


3. 7 


4 


256 


15 


656 


1312, 3936 


0.116 


64 


3, 8 


3 


192 




820 


1640, 4920 


0.093 


64 


4, 10 


3 


192 




1024 


2048. 6144 


0.074 


64 


5. 12 


2 


128 




1281 


2562, 7686 


0.059 


64 


6, 15 


2 


128 




1601 


3202, 9606 


0.048 


64 


7, 19 


1 


64 


20 


2001 


4002, 12006 


0.0381 


64 


8. 23 


1 


64 




2501 


5000, 14628 


0.031 


64 


10.28 


1 


64 



The Doppler frequency bands overlap quitre signif icantly. 
The exl:ent: of their overlap can be seen from Figure 5, 

25 where -the bands are plotted according to their respective 
velocity bins. The Doppler frequency bands are chosen to 
admit the majority of energy in a signal (>87%). The 
downsampling rates are chosen such that the time window 
over which each downsampled power estimate is formed is 

30 just greater than the duration of the typical missile 
signature* This ensures that, for each range gate 
crossing, there is a single downsampled power estimate 
which contains more than half the energy of the missile 
signature. Figure 6 illustrates this. In this figure are 

35 shown the best and worst cases, that is respectively where 
the si9:AG~L.uie, lies eri'cirfclj wiuhli. ciit "wiKit, v;indow and 
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where trhe signature spans the Junction between one window 
and the next equal ly- 

The pattern of Doppler observed over an engagement depends 
5 on the closest distance of approach of the missile to the 
receiver (the "miss distance"). This point is illustrated 
in Figure 7. This figure shows how the radial component of 
velocity varies for different scalar miss distances. The 
six different symbols represent the six range gates, the 

10 uppermost symbols on the figure representing the outermost 
range gates. In the worst case (the fourth gate out, 
represented by the full diamond symbols), the radial 
velocity can vary by a factor of 2.5 over the engagement. 
Such a change in velocity corresponds to a wandering of the 

15 signal through up to 5 velocity bins from outer gate to 
inner gate. This leads to a 4 dB signal attenuation in the 
bin detecting the missile at the outer gate compared to the 
inner gate- In the worst cases, the missile flies further 
between detections than when it has a low miss-distance and 

20 so the difference in time^^ intervals between events is 
little changed. 

Step B: Noise estimation and thresholding 

25 

At this stage in the procedure there is a set of 16 signals 
(that is, one per velocity bin) associated with each of 48 
channels ( 8 receivers with each having 6 range gates ) . As 
mentioned above under Step A, the sampling rate for each 
30 velocity bin is different. 

In Step B, an estimate is made of the background noise, a 
threshold value of noise is derived, and the signal peaks 
are compared with this threshold value. In the preferred 
35 embodiment, the signal peaks are compared with a threshold 
value which is delayed by at least two engagement times so 
that the detection of recede • events is not inhibited by the 
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higher thresholds which will be caused by the target 
activity- The relevant engagement time is the slowest 
possLhlB engagement time, in other words ^ 1.25 seconds* 

5 The easiest method of achieving a noise estimate is by what 
is termed herein the "block method". Essentially this 
involves forming a noise estimate by averaging the signal 
values over a block of samples. Each block is typically 
around 8 seconds in duration, in other words, considerably 

10 longer than twice the slowest engagement time, so that 
neither the missile approach nor recede can desensitise the 
receiving means. On the other hand, the block length is 
preferably not too much longer than * twice the slowest 
engagement time, since the expected noise is not 

15 "stationary" (in other words, it can vary quite rapidly 
with time ) . The block length is also preferably chosen so 
that in multiple missile firings a first missile does not 
so desensitise the receiving means that it prevents a 
second missile fired shortly after the first from being 

20 detected . 

When a block of signal values has been formed, the result 
is converted into a threshold by multiplying by a constant 
(say 4) and adding another small empirically derived 

25 constant (to avoid excessive sensitivity when the noise 
level is low) . The resultant threshold value is then fed 
into a delay line of a few seconds in duration, and stored 
for use at a later time. The choice of the duration of the 
delay line is made on similar principles to the length of 

30 the block ( see above ) . In order to avoid sudden steps in 
the threshold it is possible to calculate an incremental 
value which takes the threshold from the current level to 
the level to be used at the end of the next block. 

35 The above process can be implemented very efficiently 
computationally. This is because most of the work is done 
at the end of a block, say every 16 or 32 samples. 
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Normally, when eacli sample arrives, only 1:wo additions are 

required, one t:o sum -the noise power in "the current: block, 

and one to add the increment to the threshold. Provided 

the block lengths are chosen to be powers of two, any 

* 

5 divisions necessary are simple shifts. 

Because there are separate velocity bins, in principle the 
noise estimation could be optimised according to the 
velocity- In practice, because the data samples have 
10 already been downsampled in accordance with velocity, the 
optimum for each velocity is to carry out the identical 
process . 

Having defined a threshold noise value ^ each new power 
15 sample is compared with the threshold to decide whether 
there has been a detection or not at the particular range 
gate and receiver. This decision is represented., as a 
single bit of information (1 or 0). The single . bit is 
shifted into the bottom of a 16 (or possibly 32) bit word 
20 to maintain a record of the-- most recent pattern of 
detections . 

Due to the structural complexity of the above operations, 
the downsampling, noise estimation and thresholding are 
25 implemented in the preferred embodiment with a programmable 
processor. 

Step C: Processing detections from a particular receiver 
(pattern matching) 

30 

As mentioned in Section 3.3.1 above, a physically realistic 
range gate crossing sequence is hyperbolic. Such a 
sequence is shown in Figure 8. This figure provides a plot 
of the bistatic range of the missile against time- (The 
35 "bistatic range" is the sura of the ranges of the missile as 
measured between a physically separate radar transmitter 
and receiver and the missile. ) The plot is for a missile 
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travelling at a constant 180 m/s velocity and for a nearest 
miss distance of 7 metres. 

In the preferred embodiment, as mentioned above, for each 
5 of the 16 velocity bins, for each of the 6 range gates, for 
each of the 8 receivers, there is a 16 bit word defining 
the recent detection history. The next task therefore is to 
look at the recent detection history for all the range 
gates associated with a particular receiver (that is, at 

10 the evolution pattern of the detection signals) and 
determine whether the pattern of detections is indicative 
of likely engagement, thereby to determine an overall 
detection decision for the particular receiver- For 
example, in Figure 9 are illustrated two acceptable 

15 patterns attributable to an engagement that crosses all 
gates (if the gates are equally spaced). The first example 
represents a higher velocity engagement, the second 
represents a lower velocity engagement, both engagements 
being shown for the same velocity bin and therefore at the 

20 same downsampling rate. Only the approach pattern is shown. 
If a recede pattern were also included, then the pattern 
taken as an entirety would appear roughly hyperbolic. 

In order to convert the above principles into a suitable 
25 pattern matching algorithm, a set of rules has to be 
defined as to what are acceptable patterns and what are 
not. The basic rules are described below: 

1 ) The detection of interest at a particular range gate 
30 is the first one after a period of no detections. 

Subsequent detections are irrelevant, because in a 
typical engagement plenty of subsequent detections are 
to be expected, for example, due to multi-path 
reflections . 



35 



2) The detections of interest are preferably stepped by 
an amount ranging between one and two time steps per 
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range gate, -to avoid either too much or too little 
data being generated for each range gate- 



Due to the imperfections of real life, limited 
irregularities are tolerated, such as a detection at 
one range gate in advance of a correct detection 
pattern, or the absence of a detection at one range 
gate* However, the detections must in general 
progress monotonically . 

Detections at outer range gates only are regarded as 
indicative of greater miss distance. 

Note that an equivalent, time reversed, set of patterns can 
15 be defined for recede events; However, if the target 
aircraft is actually hit, then the pattern of recede 
detections is arbitrary because target fragments may be 
detected • 

20 The implementation ~ of the above ^ pattern matching,, can be 
achieved using a small amount of custom logic and a 
Programmable Read Only Memory. Alternatively, it may be 
carried out in a programmable processor. The actual 
implementation employed in the preferred embodiment is 

25 described in more detail in Section 4.6 below. 



3) 



10 



4) 



In the remainder of this sub-section, the logic necessary 
to define a gate score for each range gate (at each 
velocity bin value), and to convert the gate scores for 

30 each receiver at each velocity bin value into an overall 
detection decision and score for the individual receiver, 
for one particular velocity bin, is described. In order to 
do this, reference is first made to Figure 10, so that so- 
called "first and second points of engagement" can be 

35 defined. These points of engagement actually appear as 
lines on the detection signal pattern shown in Figure 10, 
emanating from a common point at the outermost range gate, 
and defining therefrom the slowest and fastest acceptable 
detection patterns, corresponding, respectively, on the one 
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hand, to -the slowest possible missile speed and maximiim 
miss distance and, on the other hand, to the fastest 
possible missile speed and minimum miss distance. Patterns 
falling between these extremes are deemed to be acceptable. 

The logic necessary to define a gate score for each range 
gate is as follows: 

a) If there, are any detections in the period before the 
first point of interest (that is, in the tail of the 
pattern, as illustrated in Figure 10), a gate score of 
-1 is assigned. The length of the tail needs to be 
chosen; this is typically 15 samples in the preferred 
embodiment. A negative score is given since 
detections in this region would not be consistent with 
the approach of a missile. 

b) If there ,are no detections up to and including the 
second point of engagement, the gate score is 0. 

c) If there are detections between the first and second 
points of engagement, a score is assigned which 
indicates the position of the first detection from the 
first point of engagement. 

The various possible gate scores, and the number of bits 
required to represent them, is as follows: 



10 



15 



Outer ga-te: 


-1, 


0, 


1 










2 


bits 


Gate 5: 


-1. 


0, 


1, 


2 








2 


bits 


Gate 4: 


-1, 


0, 


1, 


2, 


3 






3 


bits 


Gate 3: 


-1, 


0, 


1, 


2, 


3, 


4 




3 


bits 


Gate 2: 


-1, 


0, 


1, 


2, 


3, 


4, 


5 


3 


bits 


Inner gate: 


-1, 


0, 




2, 


3, 


4, 


5, 6 


3 


bits 



Referring to Figure 10, 



TOTAL: 16 bits 
the various gate scores basically 
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represent: the number of down- sampled periods from the 
detection at the outer gate to a detection at a closer 
gate. A score of -1 means a detection In advance of the 
line representing the first point of engagement. A score 
5 of 0 represents no detection. A score of 1 represents a 
score actually on the line representing the first point of 
engagement. A score of 2 represents a delay of one unit 
from the line representing the first point of engagement, 
a score of 3 represents a delay of two units, and so on, 

10 

The 16 bits can then be Input Into a PROM In the form of a 
look-up table or the like to determine whether the pattern 
Is acceptable or not, with the various 'gate scores acting 
as address Information for the look-up table. The output of 

15 the PROM will be "an overall detection decision score for 
the Individual receiver, for one particular velocity bin. 
(This is also termed "the instantaneous score for the 
receiver, velocity combination".) Since the PROM which is 
used has an 8 bit output, this output is used to contain 

20 both information about the quality of the engagement and 
also other miss distance information, most specifically how 
many range gates the missile has crossed. 

Ultimately, an overall assessment of the engagement is 
25 provided on the basis of the Instantaneous velocity 
combination scores for all the receivers. To allow for the 
skew between different receivers a time history of the 
scores for each receiver is stored, corresponding to the 
maximum time offset between receivers. In the present 
30 embodiment, in one down-sampled time period the missile 
would travel at least 3 m, and so this history need only be 
three samples per velocity bin per receiver. ^ 

Step D: Amalgamate results from all receivers 

35 

The Instantaneous velocity combination scores from each of 
the receivers for each of the velocity bins are amalgamated 
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t:o provide an overall assessment for the likely engagement. 

Although there are several possible ways of providing an 
overall assessment, in the preferred embodiment an 
5 engagement is scored at a particular velocity bin If the 
score, for each velocity bin, totalled across all the 
receivers is greater than six, and, to be more precise, if 
more than two receivers score more than 3 range gates or 
more than 3 receivers score two consecutive range gates in 

10 any of the sixteen velocity bins (the score being 
calculated as the number of range gates times the number of 
receivers). The latter condition (more than 3 receivers 
scoring two consecutive range gates) arises in order to 
detect a missile which only traverses the outermost range 

15 gates (say, those range gates beyond 100 feet, 30 metres). 
Approach and recede detections may be combined in this 
score but an engagement is not scored if one receiver 
scores on approach and recede only- Due allowance must be 
made for the skew between the hyperbolic range - time 

20 curves of the different receivers (the centre of the 
ellipsoids may be separated by 6 metres and so the skew of 
the curves might be this distance divided by the speed ) . 

Thus the overall engagement score for a given velocity bin 
25 is evaluated from the maximum scores from the individual 
receivers at the last three time periods in down- sampled 
time. Where an engagement is scored at multiple velocity 
bins the engagement at the lowest velocity bin is the one 
which is taken into account, since this is likely to 
30 contain the most complete information concerning the 
engagement. In practice, the lowest such velocity bin is 
usually likely to have a velocity range centred somewhat 
lower than the true missile velocity owing to the reduction 
in the missile tangential velocity as it passes the target 
35 aircraft. 

The typical result obtained from a real engagement in time 
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domain simulations is "thai: po1:ent:ial engagementrs are 
regis-tered at more than one of the assumed speeds, for 
example at the actual speed and at each of the adjacent 
speeds. These tend to be fairly consistent in the 
5 estimated time of arrival. Occasionally, apparently 
spurious engagements are registered in other velocity bins 
on signals that occur during the engagement, that is, those 
due to multipath, acoustic noise, and so on. This set of 
potential engagements can be sorted quite easily by looking 

10 at the overall engagement scores (figures of merit) for 
each one. In general, the result that yields the highest 
overall engagement score is that which most accurately 
indicates the missile's speed and time of arrival. Indeed, 
the existence of the other results can be used as 

15 supporting evidence that* a "true engagement has actually 
been detected. 

The estimate of missile speed yielded by the potential 
engagement with the highest overall engagement- score can 

20 then be confirmed or refined by comparing the- approach and 
recede results. Since the approach result yields a time of 
arrival at the outer gate, and the recede result yields a 
time of exit, another approximate estimate of missile speed 
can be made. If this is consistent with the first estimate 

25 of missile speed one can have further confidence that the 
algorithm has correctly identified the actual range gate 
crossings in both approach and recede. 

An alternative preferred embodiment takes into account more 
30 fully the fact that the Doppler varies from outer gate to 
inner gate. In this embodiment, final processing is carried 
out across some velocity bins. For example, detection 
signals from lower velocity bins at closer range gates may 
be combined with detection signals from higher bins, at 
35 farther gates. 
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4. ARCHITECTURE AND PROCESSING LOAD CALCULATIONS 

4.1 Archi-tecture 

5 The architecture of one preferred' embodiment of the ground 
station portion of the missile scoring radar system of the 
present Invention Is now described with reference to Figure 
11. This figure also shows the rate at which data Is 
transferred between the various hardware elements. 

10 

The incoming data stream received via a receiving antenna 
(not shown) on the ground station is first passed to 
Programmable Read Only Memory or programmable processor 400 
where the data which had been compressed on the target 
15 aircraft, and then transmitted down the telemetry link. Is 
decompressed. 

The decompressed data stream is passed at a data rate of 
1.68 M words per second to a 24k by 10 Buffer RAM 402, and 
20 thence, at a rate of 3.35 M words per second, to an FFT 
processor 404, described in more detail below. 

The output of the FFT processor 404 is taken at a rate of 
1.68m words per second to first and second multipliers 406 
25 and 408, whose outputs are themselves added by adder 410. 
The FFT processor 404, multipliers 406 and 408 and adder 
410 function to produce an output which represents power 
spectral densities at various frequencies. Their function 
is described in more detail in Sections 4.2 and 4.3 below. 

30 

In an alternative preferred embodiment, instead of an FFT 
processor, a processor carrying out a similar function 
(such as a wavelet transform processor) could be used. 

35 The output from the FFT processor 404 is in floating point 
format. Thus, in shifter 412 this output is scaled 
appropriately, and then buffered in a 64k by 32 Buffer RAM 
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414. The RAM 414 holds frequency bin dat:a. Data is outzput: 
from the Buffer RAM 414 at a rate of 2.4 M words per second 
to an accumulator 416, which sums the data received over a 
period of time then outputs it at a rate of 264 k words per 
5 second to an 8k by 32 Dual Port RAM 418 J which holds 
velocity bin data- The operation of these devices above is 
described in more detail in Sections 4.3 and 4.4 below. 

The RAM's 402, 414 and 418, the FFT processor 404 and the 
10 accumulator 416 are all controlled, and in some cases also 
addressed, by the state machine controller 420* This 
controller is fed with the clock output from the PROM 400. 

The output from the Dual Port RAM 418 is passed at a rate 
15 of 264 k words per second to a 68030 programmable processor 
422 from Motorola Inc., although it will be understood that 
any similar 32-bit digital processor would suffice. As 
described in. more detail in Section 4.5 below>= the 
processor 422 functions to downsample the v. data,., produce a 
20 noise estimation, and to compare the- data with - a noise 
threshold . 

The output from the processor 422 is passed at a rate of 
152 k words per second to a PROM 424 and associated logic 
25 circuitry which function as a pattern matcher in a manner 
described below in Section 4.6. 

The output from the PROM 424 is passed at a rate of 25.3 k 
words per second to a 680x0 programmable processor 426 from 
30 Motorola Inc., although again it will be understood that 
any similar processor would suffice. The processor is for 
receiver amalgamation, as described in more detail below in 
Section 4.7. 

35 In a final step in this portion of the procedure, the 
output from the processor 426, in the form of information 
on missile engagements, is passed at a low rate to a 68 k 
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processor 428* 

The PROM 400 passes calibra-tion data and t:he like "to -the 
processor 428, It: also passes, in another portion of the 
5 procedure, all the decompressed data to a 4 Mbyte Buffer 
RAM 432, which accumulates current data on the fly, 
constantly overwriting less recent data. The purpose of the 
Buffer RAM 432 is to keep a rolling memory of data so that, 
if the processor 428 recognises a potential engagement, the 
10 relevant portion of the data can be passed, under the 
control of the processor via an interface 430, to a Direct 
Memory Interface 434, thence to a Small Computer Standard 
Interface (SCSI) 436, and finally to a workstation (not 
shown), where the data can be assessed in detail. 

15 

The preferred embodiment described above would be 
relatively inexpensive to implement, but relatively 
inflexible. In another preferred embodiment, which would 
be more expensive but more flexible, ail the functions 

20 described with reference to Figure 11 are implemented on an 
array of, say, eight digital programmable signal 
processors, such as the Texas Instruments TMS 320C40, 
together with appropriate data storage facilities. In this 
event, not even the state machine controller would be 

25 required. ' 

In this case a processing engine is formed within a VME 
(trade mark) architecture comprising 8 TMS320C40 processors 
on Motherboards such as the Loughborough Sound Images (LSI) 
30 DBV42 and DBV44. These Motherboards offer 2 TMS320C40's on 
the DBV42 and up to 8 such processors on the DBV44. The 
DBV42 also has up to 16 M words of RAM for the revolving 
buffer. The processors would be allocated to functions as 
follows: - 

35 

1. Processor 1 - Data input and decompression 

2, Processors 2 to 5 - FFT or Wavelet Transform 
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3. 



Processor 6 



Velocitry bin formatrion 
Downsampling, noise 
estimatrion, thresholding and 
pat-tern matching 
Receiver analgamation 



Processor 7 



5 



5. 



Processor 8 



Various features of the preferred embodiment of the ground 
station portion of the missile scoring radar system of the 
pr&s&n-t invention are mow described in greater detail. 



4.2 FFT Processing 

15 The FFT processing is planned around a custom FFT chip, in 
the preferred embodiment the Plessey PDSP16510 stand alone 
FFT processor, as the processor 404. As discussed above, 
it could alternatively be implemented on an array of 
Digital Signal Processing chips such as TMS320C40 (Texas 

20 Instruments), or on a standard programmable processor. If 
the TMS320C40 chip were used, four processors would be 
required to implement the FFT; in this case the hardware 
would be commercial-off-the-shelf, but relatively 
expensive « 



Each group of 512 time-consecutive range gate data samples 
is processed twice, once with the FFT processor in the 2 x 
256 point mode and once with the processor in the 8 x 64 
point mode. The most efficient method of doing this is to 

30 process all channels in one mode, then process them all 
again in the other mode. This is because the FFT processor 
normally accepts input data, transforms a set of data, and 
outputs data in a parallel (pipeline) fashion. Before 
changing mode, it would be necessary to complete the 

35 retrieval of output data before commencing the input of new 
data. 



10 



still with reference to Figure 11. 



25 
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The fact that each batch of channel data needs to be 
processed twice means that it must be fed into the FFT 
processor twice. This is not a problem because it is 
necessary anyway to buffer the data before the FFT input 
5 (in the Buffer RAM 402) in order to unscramble the channels 
that arrive from the telemetry link. 

The required processing load is 48 lots of 2 x 256 point 
transforms and 48 lots of 8 x 64 point transforms, plus 

10 time to change mode twice, in 14.6 ms (where 14.6 ms 
corresponds to 512 samples at* 35000 Hz). A 2 x 256 point 
transform takes 25,8 fis in the Plessey chip, while the 8 x 
64 point transform takes 20.5 fis (assuming a clock rate of 
40 MHz). The processing time required is therefore 2.2 ms. 

15 Allowing another four periods to get data in and out at the 
changeover, this increases to 2.3 ms. Thus the chip would 
only be 16% utilised. 

4.3 Squaring of FFT Outputs 

20 

The next stage is to square and then sum the FFT outputs in 
order to convert the real and imaginary values of the 
output into power spectral densities. For this, the real 
and imaginary FFT outputs are passed to the respective 
25 multipliers 406 and 408 and thence to the adder 410 for 
squaring and adding on the fly. The outputs of this 
process are then dumped in the Buffer RAM 414 via the 
shifter 412 for subsequent resorting. The multipliers and 
adder are custom chips. 

30 

For each FFT there is one complex output for every two real 
inputs. Given that all input data is processed twice, 
there is one complex output to be processed per input 
sample. Thus, complex outputs must be processed at a rate 
35 of 48 X 35000 Hz « 1.68 MHz. Assuming that two multiplier 
chips are used, they must be able to process faster than 
590 ns each. In practice, the output data will be 
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ext:ract:ed in a more burst:y fashion and faster multiply 
rates will be required. However, since appropriate 
multipliers are readily available to less than 100 ns, for 
example a Cypress (trade mark) CY7C516 chip, there is no 
5 problem. 

4.4 Formation of Velocity Bins 

The output sequence of the frequency bins from the FFT is 
10 somewhat Jumbled and needs unscrambling. This is done as 
part of forming the velocity bins. Basically, frequency bin 
powers are read out of the Buffer RAM 414, and then summed 
together in the accumulator 416 to produce a velocity bin 
value. This value is then passed to the next process via 
15 the Dual Port RAM 418 (or else possibly via a Fast In Fast 
Out chip). The necessary addressing sequence is controlled 
by a counter and address mapping ROM ( not shown ) • 

The number of additions to be - implemented in the 
20 accumulator 416 is dependent r on the number of frequency 
bins in each velocity bin. For the bin sizes currently 
envisaged 710 adds are needed per channel per 14.6 ms time 
period, making a total rate of 48 x 710 per 14.6 ms = 2.4 
M additions per second or 420 ns per addition. This is 
25 easily achieved with a single adder. 

The number of velocity bins produced at the end of this 
process are 80 (« 2x8 (from the 2x256 point FFT processing) 
■f 8x8 (from the 8x64 point FFT processing) per channel per 
30 14.6 ms period, making a total of 264 k words per second, 
or one word per 3 . 8 //s . 

4,5 Downsampling, Noise Estimation, Thresholding 

35 Because of the requirement for downsampling, the procedure 
described in this section is structurally quite complex; it 
is therefore implemented in the preferred embodiment in a 
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programmable process . 

The "time -taken to receive 512 signal samples on each 
channel is 14.6 ms (=512/35000 ) . The following da1:a are 
5 output from the accumulator (velocity bin generator) 416: 

(a) For each of the 8 lowest velocity bins: 2 samples per 
channel =2x8x48- 768 samples, and 

10 (b) For each of the 8 upper velocity bins: 8 samples per 
channel = 8 x 8 x 48 = 3072 samples. 

The downsampling procedure is as follows*. For each velocity 
bin of each channel, the power data is downsampled (in the 
15 accumulator 416) by adding samples at successive time 
periods. As described above, the purpose of the 
downsampling is to render the data velocity independent. In 
the preferred embodiment, the downsampling factors vary 
between 2 and 8 • 

20 

When downsampling is complete, the new (downsampled) power 
estimate is compared with a noise threshold based on 
delayed low-pass filtered noise. If the power is above the 
threshold, this fact is recorded, for instance as a binary 

25 bit in a word, for later processing by the pattern matcher 
PROM 424. The noise is estimated as a long time estimate of 
received power appropriately delayed to stop the missile 
detections on approach de-sensitising the receiver. The 
noise threshold is updated periodically ( say every 8 

30 seconds ) . 

The noise estimation and noise thresholding procedure is 
described below in terms of suitable computer program steps 
(written in "pseudo-code"). 

35 

Comment: The following must happen once every 14,6 ms 
do Chan = 1, 48 
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do vbin « 1, 16 
Comment: since the "do vbln" loop Is a time-critical loop. 
It Is preferably written In In-line code (that Is, with no 
Jxmps, conditional statements or the like) 
5 do samp = 1, nsaxap(vbln) 

Comment: nsamp Is either 2 or 8 (t:orrespond±ng to the 2x256 
or the 8x64 FFT processing respectively) 
feirch ±nput:_sample 

downsum(chan, vbin) += input:_sample 

10 

increment down_counter( chan, vbin) 

if (down_count:er>=downsample_f actor (vbin) ) then 

15 Comment: a new downsampled value ls--now^ available 

Comment: update sum of powers for noise estimation 

power s downs um( Chan, vbin) 
power sum ( chan , vbin ) + » power* - 

20 

Comment: update noise tlireshold estimate 

threshold (chan, vbin) delta(chan, vbin) 

Comment: compare new value with current tlireshold 
25 Comment: "detect" Is like a delay line stored In a 16 bit 
word 

shift detect (Chan, vbin) one place left 
if (power > threshold ( chan, vbin)) then 
set bottom bit of detect(chan, vbin) 
30 endif 

output detect (chan, vbin) to pattern 
matcher 

Cozzzmentr If necessary , update noise estimation delay line 
35 increment block_count (chan, vbin) 

if(block_count(chan,vbin) >= 
block len(vbin)) then 
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Comment: sliift delay line (or change pointer to circular 
buffer) 

do i e delay_len{ vbin) , 2, -1 

delay ( Chan, vbin, i) = . delay (chan, 
5 vbln, 1-1) 

enddo 

Comment: calculate noise tJireshold based on last block and 

enter Into delay line 
10 delay(chan, vbin, 1) = 

power_suin(chan, vbin)/block_len( vbin) 
* thres_mult( vbin) + t:hresh_add( vbin) 

Comment: compute threshold Increment 
15 delta (chan^ vbin) = 

(delay(chan, vbin, delay_len{ vbin) ) - 
•threshold ( Chan, vbin) ) / 
block_len( vbin ) 

20 Comment: zero the summing variable and the counter 

power_suin(chan, vbin) = 0 
block_count:(chan, vbin) = 0 

endif 

25 down^counter = 0 

down^suin(chan, vbin) = 0 
endif 
enddo 

enddo 

30 enddo 



The above computer program steps describe the relevant 
procedure to be carried out. Of course, it is possible to 
carry the procedure out in different ways; the fastest way 
35 to implement it may be to put the channel counters as the 
innermost loops. 
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4*6 Pa-t-tern Matciiing 

Patt:ern matrchlng was described In some delrail in Section 3 
above. It Is Implemented with a small amount of custom 
5 logic (large AND gates ^ priority encoders, and so on) and 
a single PROM 424. The custom logic would probably be put 
Into a PAL or Xyllnx (trade mark) chip. 

The pattern matcher is fed with 6 words (one for each range 
10 gate) (16 or 32 bit) for each receiver and each velocity 
bin at the downsampled rate. For the fastest velocity 
bins, the downsampled rate is 35000/64 per channel, and 
thus there are 35000/64 x 48 words per second to evaluate, 
making 26,300 per second. Over all velocity bins the total 
15 is 152,000 per second.- 

Since each pattern matching operation produces two results 
per receiver (one approach, one recede), the-gvout put from 
the pattern matcher, is 15200O/6 x 2 51 k -results per 
20 second . 

Of course, it is also possible to use a more programmable 
rather than a hardware approach. In this case, a score for 
each receiver is generated from the number of consistent 
25 range gates detecting the missile. For instance, if for a 
receiver 4 range gates detect the missile consistently but 
there is one spurious event, one might generate a score 
from: - 

30 Score = 4 - IxFalse Alarm Penalty 

The scores for each receiver are then assessed by the 
receiver amalgamation process. 

35 4.7 Receiver Amalgamation 

Receiver amalgamation is. accomplished using the 
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programmable processor 426. From the above stage, new 
results are arriving at the rate of one per 20 ^s. These 
would arrive in batches of 16 related to particular 
velocity bins. The first task of the receiver amalgamation 
5 processor 426 would be to correlate the results across the 
receivers. By this means the majority of the data is soon 
discarded. In a preferred embodiment, this first step takes 
place as an interrupt driven process. A processor such as 
the 680x0 would have up to about 100 instructions per 
10 sample available to process the incoming data. 

Any remaining work involved in correlating across velocity 
bins and between approach and recede events is then carried 
. out as a background process. The process would be to 

15 identify engagements when consistent activity at a 
particular velocity bin is seen at a set of receivers. 
Typically this would occur if an engagement score threshold 
was exceeded at more than one receiver, A typical 
engagement score threshold would be 3 range gates detecting 

20 in the same velocity bin for 2 receivers. 



It will be understood that the present invention has been 
described above purely by way of example, and modifications 
of detail can be made within the scope of the invention. 
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CLAIMS 

1. Apparatus for determining a displacement characteristic 
5 of an object, comprising 

means for transmitting, at a given transmission time, 
a probe signal towards the object; 

means for receiving the probe signal returned by the 
object; 

10 means for generating a detection timing signal at a 

delay after the transmission time, corresponding to a 

selected range for the object; 

first detecting means, coupled to -the receiving means 

and .responsive to the timing signal, for detecting the 
15 returned probe signal occurring at the delay of the timing 

signal; and 

second detecting means for detecting the relative 
timing of the detection timing signal and of a^orsignal 
having a predetermined timing with respect ; to the 
20 transmission time; 

whereby the object displacement characteristic can be 
determined from the relative timing of these two signals. 

2. Apparatus according to Claim 1 wherein the transmitting 
25 means is adapted to transmit the probe signal at regularly 

repeated transmission times, at a given repetition 
frequency, and the second detecting means is adapted to 
sample the signal of predetermined timing at a sampling 
frequency lower than the repetition frequency, 

30 

3 . Apparatus according to Claim 2 wherein the sampling 
frequency is less than ICQ raillionths of the repetition 
frequency, preferably lees than 20 or even S millionths. 

35 4. Apparatus according to Claim 2 or 3 including two 
oscillators for timing the operation of the transmitting 
and second detecting means so that the repetition frequency 
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is timed by one oscilla1:or and the sampling frequency is 
timed by the beat frequency between the two oscillators. 

5. Apparatus according to Claim 1, 2, 3 or 4 wherein the 
5 transmitting means is adapted to .transmit the probe signal 

not only to the object but also via a path of fixed length, 
which does not include the object, through the ambient 
medium to the receiving means, and the second detecting 
means is adapted to detect the fixed path length signal as 
10 the signal having a predetermined timing with respect to 
the transmission time. 

6. Apparatus according to any of the preceding claims 
including splitter means for splitting the signal received 

15 by the receiving means to both the first and the second 
detecting means. 

7. Apparatus according to any of the preceding claims 
including means for • combining the signal having a 

20 predetermined timing and the detection timing signal, and 
for passing the combined signal to the second detecting 
means • 

8. Apparatus according to any of the preceding claims 
25 including means for producing a signal representative of 

the status of the apparatus in dependence on the output of 
the second detecting means. 

9 . Apparatus according to any of the preceding claims 
30 wherein the generating means is adapted to generate a 

plurality of detection timing signals, and the second 
detecting means is adapted to detect the timing of each of 
the detection timing signals relative to the signal having 
a predetermined timing. 

35 

10. Apparatus according to any of the preceding claims 
including means for determining a measure of the selected 
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range from the rela-tive "timing of -the signal having a 
predetermined timing and the detection timing signal. 

11- A method of determining a displacement characteristic 
5 of an object, comprising 

transmitting, at a given transmission time, a probe 
signal towards the object; 

receiving the probe signal returned by the object; 

generating a detection timing signal at a delay after 
10 the transmission time, corresponding to a selected range 
for the ob j ec t ; 

a first detecting step wherein the returned probe 
signal occurring at the delay of the timing signal is 
detected; and 

15 a second detecting step wherein the relative timing 

of the detection timing signal and of a signal having a 
predetermined timing with respect to the transmission time 
is detected; 

whereby the object displacement characteristic can be 
20 ' determined from -the relative^ timing of these- two^signals. 

12. A method according to Claim 11 wherein the probe 
signal is transmitted at regularly repeated transmission 
times, at a given repetition frequency, and the signal of 

25 predetermined timing lis sampled at a sampling frequency 
lower than the repetition frequency- 

13. A method according to Claim 12 wherein the sampling 
frequency is less than 100 millionths of the repetition 

30 frequency, preferably lees than 20 or even 5 millionths. 

14. A method according to any of Claims 11 to 13 wherein 
the probe signal is transmitted not only to the object but 
also via a path of fixed length, which does not include the 

35 object, through the ambient medium, and the fixed path 
length signal is detected in the second detecting step as 
the signal having a predetermined timing with respect to 
the transmission time. 
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15. A method according -to any of Claims 11 "to 14 wherein 
a plurality of detection timing signals are generated, and 
the timing of each of these signals relative to the signal 
having a predetermined timing is detected in the second 

5 detecting step* 

16. A method according to any of Claims 11 to 15 wherein a 
measure of the selected range is determined from the 
relative timing of the signal having a predetermined timing 

.10 and the detection timing signal. 

17. Apparatus for transmitting data, comprising 

input means for receiving a plurality of data 
samples ; 

15 means for evaluating a respective block data value 

characteristic of each of a plurality of data blocks, each 
data block comprising a plurality of data samples; and 

means for transmitting the respective data value for 
each block. 

20 

18. Apparatus according to Claim 17, wherein the 
evaluating means is adapted to evaluate individual data 
values representative of particular data samples in each 
block, and the transmitting means is adapted to transmit 

'25 said individual values. 

19. Apparatus according to Claim 18, wherein the 
evaluating means is adapted to evaluate the block and 
individual data values respectively as an exponent which is 

30 constant for each block and a plurality of mantissas which 
are variable within each block. 

20. Apparatus according to Claim 18 or 19, wherein the 
evaluating means is adapted to evaluate the individual data 

35 values as taking as possible values the value zero as well 
as at least one other value. 
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21, Apparatrus according to Claim 18, 19 or 20, wherein 
t:he evaluating means is adapted to evaluate the individual 
data values as taking one of an odd number of possible 
values • 

5 

22. Appara-kus according to any of Claims 18 t:o 21, 
wherein the transmitting means is adapted to transmit the 
individual data values as n bit words, each word 
representing m individual data values, the number of 

10 possible values for the individual data values taken to the 
power m being no less than 25 % lower than, preferably no 
less than 10 % lower than, and no greater than, two taken 
to the nearest power of n. 

15 23 . Apparatus according -to Claim-22 , wherein m is 3 , and 
said number of possible values is either 3 or 5* 

24. Apparatus according to any of Claims 18 to 23, 
wherein-"-^ the-^ ^evaluat ing,. means £r±3 - adapted to evaluate the 

20 individual data- value&^-in dependence ^ on tha^jOt difference 
between data samples* or samples. 

25. Apparatus according to Claim 24, wherein the 
transmitting means is adapted to . transmit at least one 

25 further data value characteristic of the absolute value of 
a given data sample or samples. 

26. Apparatus for receiving data, comprising 

means for receiving the data in the form of a 
30 plurality of block data values, each block data value being 
characteristic of a respective data block, the block 
comprising a plurality of data samples; 

means for evaluating the data samples in dependence 
on the block data values; and 
35 means for outputting the data samples. 
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27. Apparatus according to Claim 26, wherein the 
receiving means is adapted to receive individual data 
values representative of particular data samples in each 
block, and the evaluating means is adapted to evaluate the 

5 data samples in further dependence on the ihdividual data 
values . 

28. Apparatus according ' to Claim 27, wherein the 
evaluating means is adapted to evaluate the data samples in 

' 10 dependence on block and individual data values which are 
respectively an exponent which is constant for each block 
and a plurality of mantissas which are variable within each 
block. 

15 29. Apparatus according to Claim 27 or 28, wherein the 
evaluating means is adapted to evaluate the data samples in 
dependence on individual data values which take as possible 
values the value zero as well as at least one other value. 

20 30. Apparatus according to Claim 27, 28 or 29, wherein 
the evaluating means is adapted to evaluate the data 
samples in dependence on individual data values which take 
one of an odd number of possible values. 

25 31. Apparatus according to any of Claims 27 to 30, 
wherein the receiving means is adapted to receive data as 
n bit words, each word representing m individual data 
values, the number of possible values for the individual 
values taken to the power m being no less than 25 % lower 

30 than, preferably no less than 10 % lower than, and no 
greater than, two taken to the nearest power of n, and is 
further adapted to convert each received word into the m 
individual data values. 

35 32. Apparatus according to any of Claims 27 to 31, 
wherein the evaluating means is adapted to evaluate the 
data samples on the basis that the individual data values 
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have been evalualred In dependence on t:he difference be1:ween 
da-ta samples. 

33. Apparatus according to Claim 32, wherein the 
5 receiving means is adapted to receive at least one further 

data value characteristic of the absolute value of a given 
data sample or samples. 

34. Data transfer apparatus comprising apparatus 
10 according to any of Claims 17 to 25 and apparatus according 

to any of Claims 26 to 33. 

35. A method of transmitting data, comprising 
receiving a plurality of data samples; 

15 evaluating a respective block data value 

characteristic of each of a plurality of data blocks, each 
data block comprising a plurality of data samples; and 

transmitting the respective data value v for each 

block * 

20 

36. A method according to Claim 35, wherein individual 
data values representative of particular data samples in 
each block are evaluated, and said individual values are 
transmitted . 

25 

37. A method according to Claim 36, wherein the block and 
individual data values are evaluated respectively as an 
exponent which is constant for each block and a plurality 
of mantissas which are variable within each block. 

30 

38. A method according to Claim 36 or 37, wherein the 
individual data values are evaluated as taking as possible 
values the value zero as well as at least one other value. 

35 39. A method according to Claim 36, 37 or 38, wherein the 
individual data values are evaluated as taking one of an 
odd number of possible values- 
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40. A metrhod according to any of Claims 36 to 39, wherein 
the individual data values are transmitted as n bit words, 
each word representing m individual data values, the 
number of possible values for the individual data values 

5 taken to the power m being no less than 25 ^% lower than, 
preferably no less than 10 % lower than, and no greater 
than, two taken to the nearest power of n. 

41. A method according to Claim 40, wherein m is 3, and 
10 said number of possible values is either 3 or 5. , 

42. A method according to any of Claims 36 to 41, wherein 
the individual data values are evaluated in dependence on 
the difference between data samples. 

15 

43. A method of receiving data, comprising 

receiving the data in the form of a plurality of 
block data values, each block data value being 
characteristic of a respective data block, the block 
20 comprising a plurality of data samples; 

evaluating the data samples in dependence on the 
block data values ; and 

outputting the data samples. 

25 44. A method according to Claim 43, wherein individual 
data values representative of particular data samples in 
each block are received, and the data samples are evaluated 
in further dependence on the individual data values. 

30 45. A method according to Claim 44, wherein the data 
samples are evaluated in dependence on block and individual 
data values which are respectively an exponent which is 
constant for each block and a plurality of mantissas which 
are variable within each block. 

35 

46, A method according to Claim 44 or 45, wherein the 
data samples are evaluated in dependence on individual data 
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values which "take as possible values the value zero as well 
as al: least: one ether value. 

47. A method according to Claim 44, 45 or 46, wherein the 
.5 data samples are evaluated in dependence on individual data 

values which take one of an odd number of possible values. 

48. A method according to any of Claims 44 to 47, wherein 
data is received as n bit words, each word representing m 

10 individual data values, the number of possible values for 
the individual data values taken to the power m being no 
less than 25 * lower than, preferably no less than 10 % 
lower than, and no greater than, two taken to the nearest 
power of n, and each received word is converted into the m 

15 individual data'valuies. 

49- A method according to any of Claims 44 to 48, wherein 
the data samples are evaluated on the basis . that the 
individual- data values have been^^ evaluated rin dependence on 
20 the difference between data^^smiples. 

50. Apparatus for assessing the approach of an object to 
a specified location, comprising 

means for detecting the likely presence of the object 
25 in a plurality of range regions at specified ranges from 
said location, and for producing a detection signal for 
each range region in which the likely presence of the 
object is detected; 

means for determining a measure of the velocity of 
30 the object; and 

means for assessing the approach of the object in 
dependence on whether the velocity measure is consistent 
with the evolution pattern of the detection signals. 

35 51. Apparatus according to Claim 50, wherein the 
determining means is adapted to determine the velocity 
measure by sensing the velocity of the object in a 
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plurality of veloci1:y bands. 

52. Apparatus according to Claim 51, wherein the 
detecting means is adapted to produce detection signals for 

5 each velocity band, and including means for adjusting the 
evolution patterns for the velocity bands so that they are 
substantially independent of velocity. 

53. Apparatus according to Claim 51 or 52, wherein the 
10 determining means is such that the velocity bands overlap 

one another. 

54. Apparatus according to any of Claims 51 to 53, 
wherein the detecting means is adapted to detect the likely 

15 presence of the object separately at each velocity band, 
and the assessing means is adapted to assess the approach 
of the object in dependence separately on the respective 
evolution pattern of the detection signals at each velocity 
band . 

20 

55. Apparatus according to any of Claims 50 to 54, 
wherein the assessing means is adapted to assess the 
approach of the object in dependence on whether the 
evolution pattern falls between two extremes of 

25 allowability. 

56. Apparatus according to any of Claims 50 to 55, 
wherein the assessing means is adapted to assess the 
approach of the object in dependence on whether the 

30 evolution pattern approximates to a hyperbolic function. 

57. Apparatus according to any of Claims 50 to 56, 
wherein the detecting means comprxs&s a plurality of 
individual detecting means, preferably at separate 

35 locations, and the assessing means is adapted to assess the 
approach of the object in dependence on the respective 
evolution pattern of the • detection signals for each 
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individual de-tectlng means. 

58. Apparatrus according -to any of Claims 50 to 57, 
Including means for generating an estimate of background 

5 noise, and wherein the detecting means Is adapted to 
produce the detection signal according to whether the 
detecting means detects the likely presence of the object 
above a noise threshold which is dependent on the noise 
estimate . 

10 

59. Apparatus according to Claim 58, wherein the 
generating means Is adapted to generate a noise estimate 
which Is time dependent. 

15 60. Apparatus according to Claim 59, wherein the noise 
threshold is based on a delayed version of the noise 
estimate. 

61. Apparatus accorddrng .-to Claim >59raor -' 60^' wherein the 
20 noise threshold -is -basedr^onpra ..version of ^the^oise^estimate 

averaged over a given duration. 

62. Apparatus according to any of Claims 58 to 61, wherein 
the detecting means is adapted to produce the detection 

25 signal according to whether the detecting means detects the 
likely presence of the object above both of two noise 
thresholds, one of which is based on the sum of the noise 
estimate and a first constant, the other of which is based 
on a second constant times the noise estimate. 

30 

63. A method of assessing the approach of an object to a 
specified location, comprising 

detecting the likely presence of the object in a 
plurality of range regions at specified ranges from said 
35 location, and producing a detection signal for each range 
region in which the likely presence of the object is 
detected ; 
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determining a measure of the velocity of the object; 

and 

assessing the approach of the object in dependence on 
whether the velocity measure is consistent with the 
5 evolution pattern of the detection signals- 

64. A method according to Claim 63, wherein the velocity 
measure is determined by sensing the velocity of the object 
in a plurality of velocity bands. 

10 ' 

65. A method according to Claim 64, wherein detection 
signals are produced for each velocity band, and the 
evolution patterns for the velocity bands are adjusted so 
that they are substantially independent of velocity- 

15 

66. A method according to Claim 64 or 65, wherein the 
velocity bands overlap one another. 

67. A method according to any of Claims 64 to 66, wherein 
20 the likely presence of the object is detected separately at 

each velocity band, and the approach of the object is 
assessed in dependence separately on the respective 
evolution pattern of the detection isignals at each velocity 
band. 

25 

68. A method according to any of Claims 63 to 67, wherein 
the approach of the object is assessed in dependence on 
whether the evolution pattern falls between two extremes of 
allowability. 

30 

69. A method according to any of Claims 63 to 68, wherein 
the approach of the object is assessed in dependence on 
whether the evolution pattern approximates to a hyperbolic 
function. 

35 

70. A method according to any of Claims 63 to 69, wherein 
a plurality of individual detecting means are provided. 
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preferably at: separate loca-tions, and tiie, approach of -the 
object: ±s assessed In dependence on t:he respect:! ve 
evolution pattern of the detection signals for each 
Individual detecting means. 

5 

71- A method according to any of Claims 63 to 70, wherein 
an estimate of background noise Is generated, and wherein 
the detection signal is produced according to whether the 
likely presence of the object is detected above a noise 
10 threshold which is dependent on the noise estimate- 

72. ,A method according to Claim 71, wherein a noise 
estimate is generated which is time dependent. 

15 73 . A method according to Claim 72, wherein - the * noise 
threshold is based on a delayed version of the noise 
estimate . 

74. A method according to Claim 72 or 73, wherein the 
20 noise threshold is bas^ad on a version of ^the noise estimate 

averaged over a given duration. 

75. A method according to any of Claims 71 to 74, wherein 
25 the detection signal is produced according to whether the 

likely presence of the object is detected above both of two 
noise thresholds, one of which is based on the sum of the 
noise estimate and a first constant, the other of which is 
based on a second constant times the noise estimate. 
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